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ABSTRACT 


V  - 

Today's  tactical  military  air  surveillance  radars  generally  operate  in 
a  stand-alone  configuration.  The  many  performance  improvements  that  result 
when  data  from  multiple  radars  of  this  type  are  merged  have  made  such  netted 
operations  an  attractive  goal  for  many  years.  A  major  obstacle  to  achieving 
this  goal  has  traditionally  been  the  difficult)  associated  with  the  registration 
of  multisensor  data,  the  expression  of  the  data  in  a  common  coordinate  system 
free  from  errors  due  to  site  uncertainty,  antenna  orientation,  and  improper 
alignment . 

This  report  presents  the  results  of  a  modest  effort  to  develop  a  self¬ 
registration  procedure  by  which  multiple  radar  sensors  operating  in  consort 
each  calculate  the  errors  in  their  data  by  comparing  it  with  data  from  the 
remainder  of  the  system  and  then  uses  the  information  to  upgrade  performance. 

The  technique  has  been  tested  with  experimental  data  and  appears  quite  capable 
of  improving  system  performance,  measured  in  terms  of  residual  inter-site  bias 
errors,  by  almost  a  factor  of  one  hundred. 
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INTRODUCTION 


I  . 

Since  an  air  surveillance  radar  was  first  developed  prior  to  World  War  II, 
military  planners  have  given  much  attention  to  the  netting  of  individual  sensors; 
the  combining  of  data  from  multiple  sensors  to  provide  a  universal  picture  of 
the  overall  air  situation.  Attempts  at  such  integration  have  traditionally 
met  with  limited  success. 

There  are  many  reasons  why  such  integration  is  desirable,  indeed  necessary 
in  some  situations.  As  weapon  system  technology  advances  and  reaction  tines 
diminish,  many  of  these  are  becoming  critical. 

A .  Quality  of  Coverage 

Individual  radar  sensors  are  limited,  both  in  the  volume  of  airspace 
(especially  at  low  altitudes)  they  can  search  and  in  the  quality  of  the  air 
track  data  they  can  form,  particularly  with  respect  to  tracking  maneuvering 
aircraft  through  rapidly  varying  flight  profiles.  Effective  integration  of 
multisensor  data  could  fill  coverage  gaps  and  could  more  rapidly  discriminate 
between  maneuvers  and  surveillance  eiiurs. 

B.  Reliability 

When  a  stand-alone  radar  site  is  attacked  successfully  all  users  of  its 
output  data  are  effectively  denied  information.  An  integrated  system  whose 
output,  is  shared  among  all  users  would  produce  data  whose  quality  would  degrade 
only  gradually  as  individual  sites  were  eliminated. 

C .  Coord i nation  of  Identity  Data 

An  integrated  air  surveillance  system  would  necessarily  employ  a  common 
track  data  base  to  which  additional  information,  particularly  information  re¬ 
garding  friend  or  foe  identity,  could  be  readily  added  no  matter  how  derived. 

A  single  positive  identification  of  a  particular  target  could  serve  to  provide 
all  data  users  with  its  identity'  status  for  the  dur.  tion  of  its  lifetime  within 
the  system,  thus  reducing  the  need  for  a  highly  reliable,  new,  IFF  system. 
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D. 


Increased  Reliance  on  Automation 


The  quality  of  sensor  data  available  from  today's  radar  equipment,  especially 
in  liCM  situations,  is  such  that  automated  processing  of  that  data  frequently 
yields  only  marginal  performance  and  human  intervention  is  necessary  to  determine 
which  tracks  are  valid,  which  targets  represent  actual  aircraft,  and  so  forth. 

Ior  man)'  years  attempts  have  been  made  to  develop  systems  that  are  less  man¬ 
power  intensive,  both  by  improving  radar  data  quality  and  by  improving  the 
sophistication  of  the  processing  algorithms.  The  capability  to  operate  on 
data  derived  from  multiple  sites  could  lead  to  significant  improvements  in 
automated  system  performance  and  result  in  a  far  higher  degree  of  automation 
using  today's  radar  systems  than  can  be  obtained  when  those  systems  are  operated 
in  the  present  stand-alone  mode. 

E.  ECM  Resistance 

An  integrated  system  would  be  naturally  resistant  to  ECM  since  the  sectors 
of  individual  radar  sensors  that  are  disrupted  by  jamming  would  effectively 
be  filled  in  by  other  radars  whose  geometry,  relative  to  the  ECM  source,  would 
be  different.  In  addition,  the  communication  links  and  the  ability  to  operate 
from  a  common  data  base  that  form  the  basis  of  the  integrated  system  would 
readily  support  a  unified  approach  to  emission  control  and  other  ECCM  techniques. 

F •  Integration  with  Other  Systems 

Other  systems  which  establish  or  employ  position/locaticn  information 
(e.g.,  JilbS,  PLRS,  inertial  navigators)  could  readily  be  integrated  with 
the  netted  radar  surveillance  system  since  its  common-grid  coordinate  system 
and  the  registration  process  essential  to  the  establishment  of  that  system 
would  readily  support  such  integration. 

These  attributes  have  been  recognized  as  advantageous  in  many  applications 
for  several  years,  and  numerous  attempts  have  been  made  to  achieve  them  by  in¬ 
tegration  of  multiple  radar  sensors.  Some  of  these  have  been  successful,  some 
have  not,  depending  mainly  on  the  position  accuracy  required  from  the  integrated 
system.  In  general,  integrated  systems  in  which  required  accuracy  is  less  than 
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a  few  miles  have  not  been  readily  achievable.  Studies  of  past  history  in  this 
area  suggest  that  the  reason  for  this  is  inadequate  ability  to  perform  registra¬ 
tion  on  the  multiple  input  sources;  that  is,  to  align  them  to  a  common  coordinate 
system.  Performance  in  this  area  has  been  traditionally  limited  by  data  process¬ 
ing  capability,  the  several  coordinate  conversion  "shortcuts"  necessary  to  per¬ 
form  processing  within  the  capabilities  of  available  digital  processors  resulted 
in  worst-case  accuracies  of  several  miles.  With  the  advent  of  powerful  and  in¬ 
expensive  real-time  computers  it  is  appropriate  to  reexamine  the  registration 
process  to  determine  whether  it  can  result  in  significant  improvements  in  accuracy. 
This  report  does  that  and  concludes  that  accuracy  improvements  of  more  than  two 
orders  of  magnitude  are  realizable  with  presently  available  technology.  A  highly 
precise  registration  algoritlun  is  .  veloped  an  l  verified  with  experimental  data. 
This  algorithm  could  form  the  basi  r. '  an  automated  self -l  egistrat  ion  system  that 
would  make  feasible  the  highly  accut  .c  integrated  air  surveillance  system  of 
the  future. 

II.  BACKGROUND 

The  earliest  air  search  radar  integration  process  was  developed  more  than 
40  years  ago  around  the  "Home  Chain"  radar  sensors  wnich  alerted  British  defenses 
to  German  air  attack  during  the  Battle  of  Britain.  System  accuracies  were  poor 
(on  the  order  of  several  miles) ,  areas  of  coverage  overlap  were  few  (outnumbered 
by  areas  of  no  coverage) ,  and  processing/data  integration  functions  were  es¬ 
sentially  manual.  On  the  other  hand,  leaction  time  and  vectoring  accuracy  re¬ 
quirements  were  sufficiently  low  so  the  system  proved  quite  effective. 

The  development  of  SAGE  in  the  years  following  World  War  II  provided  the 
first  opportunity  to  merge  automatic  computing  technology  with  radar  data  and 
brought  about  the  first  attempts  to  systematically  define  and  overcome  registra¬ 
tion  problems.  The  procedures  that  were  developed  in  SAGE  for  dealing  with 
data  on  a  single  target  gathered  by  multiple  radars  are  essentially  those  found 
today  in  many  applications. 
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For  each  volume  of  airspace  a  particular  sensor  is  defined 
as  "primary".  Other  sensors  are  defined  as  "secondary",  "tertiary", 
etc . 

An>'  track  within  that  volume  is  updated  exclusively  from  data 
obtained  by  the  primary  sensors.  If  the  "quality"  of  that  data 
(measured  in  terms  of  the  blip-scan  ratio  or  detection  probability) 
falls  below  a  threshold,  input  is  switched  to  the  secondary  sensor, 
etc . 

At  boundaries  between  volumes  for  which  different  sensors  provide 
primary  coverage,  large  discontinuities  in  tracks  occur.  Tracking 
software  is  programmed  to  ignore  these,  and  operators  are  trained  to 
live  with  them. 

This  procedure  is  employed  in  today's  version  of  SAGE,  and  is  essentially 
that  used  in  the  FAA's  enroute  radar  tracking  system,  NAS  Stage  A.  (It  should 
be  noted  that  in  the  FAA's  application  the  procedure  more  than  suffices  since 
aircraft  1 ly  prescribed  tracks  and  the  measured  parameter  of  most  interest  is 
position  along  track  versus  time  as  opposed  to  the  position  of  the  track  itself.) 

The  SAGE  procedure  for  multisensor  registration  was  adopted  by  the  Navy  for 
use  in  its  NTDS  (Naval  lactical  Data  System)  which  is  an  automated  system  for 
the  netting  of  surveillance  data  and  dissemination  of  tactical  command  and 
control  information  employed  on  most  larger  Naval  combatant  vessels.  The  sys¬ 
tem,  of  course,  involves  sensors  whose  relative  positions  vary  with  time,  and  a 
more  sophisticated  registration  process  is  in  order.  Over  the  20  years  in  wh’ch 
NTDS  has  been  operational  the  registration  software  has  improved  considerably 
through  a  sequence  of  "gr i d- lock ing”  algorithms  and  associated  algorithms  for 
the  alignment  of  data  from  multiple  radars  on  t lie  same  platforms.  Overall 
registration  accuracy  is  still  only  on  the  order  of  a  mile  or  so,  but  this  would 
appear  sufficient  for  the  purposes  of  the  system  as  it  is  currently  employed. 

Many  other  attempts  at  netting  radar  data  have  been  made  in  the  context  of 
recent,  specific  systems  such  as  the  Army's  Missile  flinder  (TSQ-73) ,  the  Air 
Force's  Tactical  Air  Control  System,  and  the  Marine  Air  Command  and  Control 
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System.  These  systems,  for  t he  most  part,  Jo  not  attempt  to  capital  ire  on  the 
advantages  inherent  in  multi  sensor  operation;  the  limited  data  integration  of 
which  they  arc  capable  is  not  considered  of  paramount  importance  and,  hence,  the 
registration  approaches  they  employ  are  unsophisticated,  in  some  instances  being 
merely  19S0  equivalents  of  the  SAG!:  approach  to  registration.  Kith  the  real  1  ra¬ 
tion  of  the  advantages  to  be  gained  by  a  more  systematic  and  comprehensive 
approach  to  data  integration,  and  the  recognition  of  the  vulnerability  of  es¬ 
sentially  stand-alone  systems,  this  situation  is  slowly  changing  but  no  examples 
of  improved  registration  capabilities  consistent  with  that  change  are  presently 
apparent  in  the  field. 

Ill.  TUI:  INTEGRATED  AIK  SUKVf. Il.LANCh  SVSThM  Of  THi:  FIJI  UKP 

Present  air  search/surveillance  resources  are  numerous  and  varied  in  capa¬ 
bility.  Rather  than  considering  the  development  of  a  new  integrated  air  sur¬ 
veillance  system  using  a  new  fami.H  of  radar  sensors,  it  appears  economically 
attractive  to  consider  the  integration  of  sensors  that  are  presently  in  use,  for 
the  most  part  in  stand-alone  operation.  These  sensors  are,  for  the  most  part, 
manual  or  semi -automated ;  the  quality  of  the  data  they  produce  is  generally 
insufficient  to  support  fully  automated  track  initiation  and  updating.  It  is 
anticipated  that  the  proper  form  of  netting,  allowing  the  sharing  of  raw  target 
data,  will  support  a  higher  level  of  automation  than  presently  realised.  Some 
of  the  most  recent  radar  equipments,  employing  fairly  sophisticated  radar  signal 
processing  techniques,  provide  data  of  sufficient  quality  to  support  this  higher 
degree  of  automation.  The  trade-offs  between  increased  netting  and  this  in¬ 
creased  processing  sophistication  as  means  to  accomplishing  full)'  automated 
tracking  are  not  well  understood,  but  it  would  seem  that  material  simplification 
of  the  basic  radar  systems  could  result  from  proper  application  of  netting  and 
this  would  offset  some  of  the  expenses  associated  with  netted  operation. 

for  a  variety  of  reasons,  particularly  for  minimum  cost  and  to  minimize 
vulnerability  to  physical  attack,  the  integrated  system  should  be  distributed; 
data  merging  and  processing  should  be  performed  at  many  sites  rather  than  con¬ 
centrated  into  a  hierarchical  structure.  It  would  be  appropriate  to  make  those 
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sites  the  radars  themselves.  Suitable  addition  of  digital  processing  and  inter 
communication  equipments  to  the  radars  would  suffice  to  transform  present-day 
sensors  into  the  integrated  future  system.  Consistent  with  this  system  archi¬ 
tecture,  individual  sensor  sites  would  "broadcast"  data  on  targets  detected  and 
tracks  in  process,  and  would  employ  data  received  over  "broadcasts"  from  these 
sensors  as  well  as  their  own  detection  data  as  input  to  their  tracking  process. 

The  definition  and  development  of  the  appropriate  data  merging  and  tracking 
processes  represent  a  major  task;  it  is  clear  that  proper  data  registration  is 
an  essential  prerequisite  to  these  processes.  The  registration  process  must 
transform  position  data  on  all  targets  and  tracks  to  a  form  suitable  for  input 
to  each  tracking  process.  To  be  consistent  with  the  tactical  situation,  regis¬ 
tration  should  be  a  continuing  process.  Upon  initial  setup  of  a  network  of 
radar  sites  the  registration  process  should  automatically  begin  to  tie  the  data 
from  those  sites  together  by  appropriate  adjustment  of  site  parameters  based 
entirely  on  data  derived  from  targets  of  opportunity.  As  more  targets  with 
different  geometries  become  processed  by  the  system,  the  registration  process 
should  continually  refine  its  estimates  of  the  biases  associated  with  the  in¬ 
dividual  sites  and  their  error  mechanisms.  In  the  case  of  mobile  or  readily 
t ransportable  sites  the  process  should  support  the  tracking  of  the  position  of 
the  individual  sites  and  maintain  registration  throughout  movement  and  recon¬ 
figuration  activities. 

The  registration  issues  addressed  in  this  report  represent  the  first  step 
in  the  development  of  such  an  automated  self-registration  process;  at  the  heart 
of  any  such  process  is  an  understanding  of  the  mechanism  by  which  various  error 
sources  contribute  to  total  registration  error.  Section  IV  covers  these  in 
deta  i  1 . 

Given  a  thorough  enumeration  of  registration  error  sources,  a  self- 
registration  algorithm  to  characterize  those  sources  and  correct  data  for  their 
effects  can  be  developed;  this  is  done  in  Section  V.  Some  experimental  verifi¬ 
cation  of  the  error  model  and  se 1 f -regi st rat  ion  algorithm  are  presented  in 
Section  VI.  Section  VII  presents  conclusions  regarding  the  pertinence  of  this 
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ontiic  topic  to  an  automated  integrated  surveillance'  system  and  recommend:, 
additional  work  in  the  tracking  and  surveillance  areas  needed  to  bring  such  a 
system  into  being. 

IV.  bQi IKCI.S  01  RlulSTRAl  ION  1.KK0K 

When  a  single  air  surveillance  sensor  is  employed  for  air  traffic  control 
purposes  only,  relative  aircraft  positions  and  courses  arc  important.  Offsets 
in  range  and  azimuth  which  apply  equal 1>  to  all  aircraft  do  not  affect  per¬ 
formance.  When  it  is  desirable  to  combine  two  or  more  overlapping  sensor  out¬ 
puts,  the  situation  changes  and  a  variety  of  error  sources  must  be  considered 
(see  Table  I).  All  error  sources  which  would  contribute  more  than  a  few  meters' 
error  in  the  registration  of  aircraft  reports  between  two  radars  are  listed  in 
the  tabic  and  discussed  below. 

A .  Range  Lrrors 

There  arc  four  types  of  range  errors  due  to  range  offset,  range  clock  rate 
error,  propagation,  and  erroneous  slant  range  correction.  Range  offset  refers 
to  a  common  increment  in  range  added  to  all  range  measurements.  The  digital 
range  counters  found  in  modern  air  surveillance  radars  must  be  properly  zeroed 
to  eliminate  range  offset  by  using  an  accurately  surveyed  radar  target.  Once 
adjusted  there  is  littie  likelihood  of  further  adjustment  even  when  the  radar 
is  moved.  Errors  in  the  range  clock  can  produce  errors  proportional  to  range. 
Another  range  correction  is  required  due  to  the  presence  of  the  troposphere 
which  influences  the  velocity  at  which  the  radar  signals  propagate.  This 
error  is  a  nonlinear  function  of  aircraft  range  and  height.  Appendix  A  contains 
a  table  suitable  for  use  with  standard  atmospheres  and  radars  located  near  sea 
level.  From  the  chart  wc  see  that  corrections  up  to  over  100  meters  are  re¬ 
quired  for  long-range  aircraft. 

Under  range  error  sources  we  have  also  listed  the  obvious  error  if  slant 
range,  as  measured  by  the  sensor,  is  used  as  a  horizontal  range.  A  large  error 
is  produced.  For  instance,  at  30°  elevation  angle  the  error  in  horizontal  range 
is  about  one-quarter  of  the  aircraft's  height  (about  14  percent  of  its  range). 
Accurate  registration  must  properly  account  for  aircraft  height.  How  this  is 
accomplished  is  explained  below  under  coordinate  conversions. 


SENSOR  REGISTRATION,  SOURCES  OF  ERROR 


COORDINATE  CONVERSION:  ACCURATE  ALGORITHMS 


B.  Azimuth  brrors 


Trrors  in  azimuth  can  be  caused  by  incorrect  alignment  of  the  radar 
antenna  with  the  reference  azimuth  (i.e.,  north)  and  by  incorrect  alignment  of 
its  axis  of  rotation  with  the  local  vertical.  Azimuth  offset  error  occurs  when 
the  radar  does  not  point  to  true  north  when  its  azimuth  readout  indicates  north. 
North-seeking  gyro  compasses  are  available  which,  when  left  stationary,  will 
produce  a  true  north  reading  to  2  milliradians  within  about  four  minutes.  These 
would  be  accurate  enough  for  use  with  most  mobile  radars.  Alternatively,  initial 
alignment  using  a  magnetic  compass  could  be  followed  by  fine  alignment  using  the 
self-registration  algorithm  described  later. 

Another  possible  source  of  azimuth  error  is  antenna  tilt.  Tor  an  ordinary 
2D  (i.e.,  non-height  finding)  rotating  sensor  with  horizontal  azimuth  bearing 
two  types  of  antenna  tile  can  occur.  The  entire  azimuth  bearing  (and  antenna) 
may  be  tilted  or  the  azimuth  bearing  may  be  level  but  the  electrical  axis  of  the 
antenna  may  be  tilted.  Both  tilts  produce  errors  proportional  to  the  tangent  of 
the  elevation  angle  of  the  target.  Bearing  tilt  will  produce  errer  terms  pro¬ 
portional  to  the  sine  and  cosine  ot  the  azimuth  angle  (see  error  equations  in 
Section  V).  Azimuth  bearing  tilt  can  be  minimized  either  by  carefully  leveling 
the  azimuth  bearing  or  by  using  tiltmeter  outputs  to  correct  aircraft  position 
reports. 

Tilt  of  the  electrical  a/.is  above  the  azimuth  bearing  should  be  carefully 
calibrated  at  the  factory  and  either  removed  or  used  to  correct  airciaft  posi¬ 
tion  reports.  Steerable  3D  sensors  may  possess  other  tilt  errors  depending  on 
the  arrangement  of  their  azimuth,  elevation,  and  traverse  axes. 

C.  Time  Errors 

Two  possible  sources  of  registration  error  are  connected  with  the 
timing  of  aircraft  position  reports.  Various  sensors  may  be  offset  in  time 
between  one  another.  Also,  some  sensors  may  have  clocks  that  run  fast  or  slow. 

A  method  for  accurately  setting  the  sensor's  clock  should  be  devised  using  the 
intersensor  communications  network  and  accurate  crystals  should  be  used  in  all 
clocks.  Ter  registration  accuracy  of  10  meters  on  Mach-2  targets  the  clocks 
should  be  accurate  to  20  msec.  Present  practice  is  to  assume  that  the  time  of 
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the  report  over  the  communications  link  is  the  time  at  which  the  aircraft  posi¬ 
tion  was  observed.  Since  more  than  a  20-msec  variation  may  exist  in  reporting, 
a  time  tag  should  be  incorporated  within  every  position  report. 

D .  Radar  Location  Errors 

Registration  errors  will  also  exist  if  sensor  locations  are  in  error. 

A  number  of  navigation  systems  can  be  used  to  provide  accurate  location  data. 
These  are:  the  Position  Location  and  Reporting  System  (PLRS) .  the  Joint  Tactics’ 
Information  Distribution  Systems  (JTIDS) ,  the  Global  Positioning  System  (GPS), 
and  various  inertial  navigation  systems.  We  will  not  review  each  of  these 
navigation  systems.  All  but  the  last  claim  to  provide  location  accuracies 
comparable  to  the  desired  10  meters. 

E .  Coordinate  Conversion  Errors 

Errors  can  also  result  from  inaccuracies  in  the  process  of  converting 
target  coordinates  to  system  coordinates  if  traditional  approximations  are  made 
and  all  corrections  are  not  included.  In  this  regard,  choice  of  the  proper 
coordinate  system  is  important. 

Sc-veiai  good  reasons  can  be  found  for  use  of  an  accurate  earth-referenced 
coordinate  system.  Such  a  system  is  also  desirable  when  other  locating  systems 
are  employed  to  locate  certain  targets  or  components  of  the  air  surveillance 
sensors  to  serve  users  also  located  in  geodetic  coordinates.  For  instance,  most 
military  aircraft  carry  inertial  navigation  systems  (INS).  When  accurately 
registered  in  geodetic  coordinates  the  surveillance  data  could  be  relayed  to 
friendly  aircraft  to  update  their  INS.  Alternately,  the  friendiy  aircraft  could 
report  their  INS-derived  positions  for  correlation  with  surveillance  sensor  data 
leading  to  positive  identification. 

When  air  defense  batteries  are  located  in  geodetic  coordinates,  it  is 
most  convenient  to  supply  target  coordinates  in  geodetic  coordinates. 

Registration  in  geodetic  coordinates  allows  easy  reference  to  a  wider 
variety  of  maps  and  accurate  placement  of  map  data  on  system  displays.  Airborne- 
radars  and  direction  finders  utilize  INS  for  short-term  position  location.  In 
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order  to  enter  their  detection  and  track  data,  the  system  should  use  geodetic 
coordinates . 

Future  navigation  systems  such  as  GPS  and  jTIDS  expect  accuracies  cn  the 
order  of  10  to  30  meters.  The  coordinate  system  and  conversion  algorithms  used 
in  an  integrated  air  surveillance  system  should  match  these  accuracy  numbers. 

This  rules  out  most  of  the  techniques  used  for  coordinate  conversion  in  present- 
day,  integrated  air  surveillance  systems. 

Present-day  systems  employ  stereographic  projection  (refs.  1-4).  The 
sensor  reports  are  projected  onto  a  plane  local  to  the  radar  and  tangent  to  a 
sphere.  When  received  by  the  master  site  they  are  transformed  onto  the  master 
plane  rangent  to  a  sphere  at  the  master  site.  A  second-order  transformation 
is  usually  employed  which  is  accurate  to  a  few  meters  (ref.  4)  over  several 
hundred  miles.  Inaccuracies  of  a  kilometer  or  so  are  typically  introduced  how¬ 
ever  in  the  stercographic  projection  to  the  local  plane.  Highly  accurate, 
direct  conversion  of  sensor  data  to  geodetic  coordinates  is  fairly  simple  (sec 
Appendix  B)  and  avoids  the  necessity  for  further  transformation  to  each  master 
plane . 

Tne  simplicity  of  data  sorting  and  the  absence  of  any  required  coordinate 
conversion  are  strong  factors  favoring  the  reporting  of  data  in  geodetic  co¬ 
ordinates.  Lvery  source  or  user  of  data  can  act  independently  in  receiving  and 
sorting  out  data  which  may  be  useful  to  him  from  data  being  broadcast  by  other 
data  sources.  Thus,  a  system  of  sensors  and  users  can  grow  to  as  large  a  net¬ 
work  as  desired  or  can  contract  to  a  set  of  autonomous  sensors.  The  use  of 
geodetic  cooidinates  provides  a  common,  well  understood  language  allowing  easy 
communication  amongst  diverse  data  sources  and  users. 

The  considerations  listed  above  suggest  that  geodetic  coordinates  are  the 
best  choice  in  a  system  of  netted  air  surveillance  sensors.  Processing  com¬ 
plexity  has  limited  their  use  in  the  past,  but  with  today's  computer  capabilities, 
their  innate  accuracy  justifies  the  slight  additional  workload  required  with  a 
geodetic  system. 

Since  situation  displays  are  generally  flat  and  the  world  round  (nearly), 
some  type  of  projection  must  be  employed.  Where  the  displayed  area  is  no  larger 
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than  one  or  two  hundred  miles  in  extent  the  Mercator  projection  is  acceptable. 

A  Mercator  projection  is  used  for  all  of  the  situation  figures  in  this  report. 
Aircraft  positions  in  geodetic  coordinates  are  easily  converted  to  a  Mercator 
projection  since  lines  of  constant  latitude  and  longitude  are  rectangular 
straight  lines  and  bearing  angles  are  accurately  preserved.  For  displays  of  very 
large  areas  some  other  projection  such  as  the  Lambert  conical  projection  may  be 
desired.  To  preserve  accuracy  and  maintain  maximum  convenience  all  aircraft 
positions  and  track  projections  should  be  calculated  in  geodetic  coordinates. 

V .  SELF-REGISTRATION  ALGORITHM 

Each  of  the  various  registration  error  sources  discussed  above  can  be 
minimized  by  more  careful  Siting,  calibration,  etc.  An  alternative  which 
would  appear  more  attractive  in  a  tactical  situation  would  be  a  computer  algorithm 
which  automatically  brings  each  sensor  into  registration  with  the  remainder  of 
the  network  whenever  its  situation  is  changed.  In  any  event,  an  algorithm  is 
required  to  assess  the  sensor's  alignment  with  respect  to  the  system. 

Upon  initial  set-up  the  best  estimates  of  position,  north  reference,  etc., 
would  be  fed  to  the  algorithm  which  would  then  examine  target  reports,  compare 
them  with  those  obtained  from  the  remainder  of  the  system,  and  automatically 
vary  the  various  registration  error  correction  values  to  bring  the  radar  into 
alignment  with  the  rest  of  the  system.  The  process  would  use  targets  of  op¬ 
portunity  and  would  presumably  converge  on  the  optimal  set  of  bias  error  cor¬ 
rections  as  the  number  of  target  detections  employed  increased  and  as  the 
entire  variety  of  target/sensor  geometries  was  satisfied.  This  section  dis¬ 
cusses  such  an  algorithm,  developed  at  Lincoln  Laboratory,  which  solves  auto¬ 
matically  for  antenna  pointing  misa 1 igament ,  range  offset  and  range  clock  error, 
time-base  differences,  and  sensor  location  errors. 

While  ideally  such  an  algorithm  would  operate  to  lock  one  sensor  into  an 
established  system  of  sensors,  the  one  discussed  here  functions  to  register  two 
sensors  with  one  another. 

!n  discussing  tne  algorithm  it  is  assumed  that  the  observations  have  been 
made  by  two  sensors  with  overlapping  coverage  for  an  extended  period  of  time  and 
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have  undergone  all  necessary  preprocessing  as  described  in  Appendix  C,  and  that 
centroided  target  reports  (range,  azimuth,  altitude,  and  time)  are  available 
from  both  sensors  for  each  aircraft  under  consideration.  The  steps  performed  by 
the  self-registration  algorithm  are: 

(1)  smoothing  and  coordinating  the  timing  of  the  target 
reports  from  each  sensor  to  generate  a  smaller  set 
of  reports  (referred  to  as  ''superpoints") , 

(2)  adjusting  the  superpoints  to  account  for  known 
biases , 

(3)  correcting  the  superpoint  ranges  for  atmospheric 
refraction, 

(4)  computing  the  covariance  matrix  cf  the  difference 
of  positions  as  observed  by  the  two  radars  for 
each  superpoint, 

(5)  solving  for  the  set  of  bias  components  using  a 
maximum  likelihood  approach,  and 

(6)  estimating  the  variances  of  the  biases. 

Step  (5),  being  iterative,  is  the  least  straightforward  and  most  computa¬ 
tionally  demanding  of  the  six.  The  steps  are  discussed  in  detail  in  the 
following  paragraphs. 

Step  1 : _ S  nothing 

Proper  bias  estimation  using  data  from  two  sensors  observing  the  same 
aircraft  necessitates  target  reports  which  are  coincident  in  time  according  to 
the  clocks  at  each  sensor  and  which  are  also  free  from  non- systematic  errors 
such  as  might  arise  from  noisy  or  garbled  target  reports.  Since  time  errors 
have  been  included  in  the  bias  model,  it  is  necessary  also  to  have  at  least 
first-order  approximations  to  the  time  derivatives  of  Tange  and  azimuth. 

Curve  fitting  (ref,  10)  with  fv  rth-degree  orthogonal  polynomials  was  found 
to  satisfy  the  requirements  listed  above.  Short  track  segments,  lasting  about 
one  minute,  where  chosen  for  fitting  to  ensure  that  aircraft  maneuvers  could  be 
adequately  followed. 
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A  computer  program  and  several  associated  subroutines  were  written  to 
accomplish  the  following  individually  for  range,  azimuth,  and  altitude: 

(a)  Identify  track  segments  for  curve  fitting  only  if  the 
associated  time  intervals  for  the  two  sensors  are  co¬ 
incident  for  at  least  a  specified  number  of  seconds 
(typically  between  10-12  observations). 

(b)  Compare  the  mean  squared  value  of  the  residuals  after 
fitting  with  a  threshold  defined  for  each  component  based 
on  expected  random  measurement  errors.  If  the  mean 
squared  value  exceeds  the  threshold,  then  the  observation 
having  the  largest  deviation  from  the  least  squares  poly¬ 
nomial  is  discarded  and  the  remaining  data  is  fitted 
again.  Should  the  mean  squared  value  still  be  too  high, 
the  time  interval  is  skipped  over  for  both  sensors  and 
control  is  passed  back  to  step  (a). 

(c)  Evaluate  the  polynomials  at  a  common  mid-point  time.  The 
resulting  "superpoint"  includes  time,  range,  range  rate, 
azimuth,  azimuth  rate,  and  altitude. 

Step  2:  Initial  Correction  for  Assumed  Biases 

In  some  cases  approximate  bias  values  may  be  known  beforehand.  For  ex¬ 
ample,  by  plotting  altitude  during  times  when  the  aircraft  is  performing  a  steep 
ascent  maneuver,  it  may  be  possible  to  estimate  a  time  offset  between  two  sen¬ 
sors  (refer  to  Figure  C-2). 

The  formulas  used  to  adjust  the  superpoints  for  initial,  assumed  bias  values 
and  also  for  the  iterative  bias  estimation  procedure  in  Step  5  are  written  below. 
SENSOR  1 

P1  =  (p  ♦  ACTD)  (1  RC1 ) 

=  0:  +  AZll  tan  •  (AZ21  ♦  AZ31  sin@1  +  AZ41  cosOj) 

SENSOR  2 

At  =  (t  +  TB1 2)  (1  ♦  TB22)  -  t 
p  2  =  (p  +  R2  +  ACTD  +  p*2  At)  (1  +  RC2) 

%2  =  e2  ♦  AZ12  ♦  e2  At  ♦  tan  n2  •  (AZ22  +  AZ32  sin  fc>2  ♦  AZ42  cos  6  ) 

^2  -  $2  +  LAT2 

^2  -  A 2  ♦  LON2 
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where:  the  subscripts  1  and  2  refer  to  sensors  1  and  2, 

Pj  and  p?  are  measured  aircraft  ranges, 

0^  and  0^  are  measured  aircraft  azimuths, 

and  r ^  are  measured  aircraft  elevation  angles, 

4<2  and  are  the  latitude  and  longitude  of  sensor  number  2,  and 

the  tilde  above  certain  quantities  identify  the  quantity  after 
correction  using  the  biases,  while  a  dot  implies  a  time  derivative. 

The  biases  in  the  above  equations  are  defined  as  follows: 

ACTD  is  the  sum  of  the  aircraft's  transponder  time  delay  measured  in 
range  units  plus  any  delay  in  sensor  1.  It  is  a  constant  for  any  one  aircraft, 
but  varies  from  aircraft  to  aircraft.  For  skin-painting  radars  it  is  a  constant. 

RC1  and  RC2  are  range  rate  errors  caused  by  inaccurate  range  clock  rates, 

AZ11  and  AZ12  are  azimuth  offset  biases, 

A221  and  AZ22  are  tilt  biases  above  the  azimuth  bearing, 

AZ31 ,  AZ4 1 ,  AZ32  and  AZ42  are  biases  describing  the  tilt  of  the  azimuth 
bearing, 

TB12  is  the  time  offset  bias  of  sensor  2  relative  to  sensor  1, 

TB22  is  the  clock  rate  difference  of  sensor  2  relative  to  sensor  1, 

R2  is  the  range  offset  bias  of  sensor  2  relative  to  sensor  1,  and 

LAT2  and  L0N2  are  the  errors  in  the  assumed  location  of  sensor  2. 

Step  5:  Atmospheric  Refraction  Correction 

Because  of  delay  caused  by  the  earth's  atmosphere  the  observed  target  range 
will  appear  to  be  greater  than  it  actually  is.  The  refraction  correction  employed 
is  based  on  a  20  x  20-element  lookup  table  characterizing  an  idealized  standard 
atmosphere  (see  Appendix  A).  Logarithmic  interpolation  has  been  used  for  values 
not  directly  found  in  the  table. 

By  way  of  example,  the  range  correction  for  a  target  with  an  elevation  angle 
of  5°  at  a  range  of  100  km  is  approximately  108  jr.  A  few  tens  of  meters  is  more 
typical  of  corrections  to  target  ranges  encountered  in  this  study. 
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For  the  experimental  data  described  in  Section  VI  no  attempt  was 
obtain  a  refraction  correction  table  corresponding  to  the  atmospheric 
the  actual  sensor  locations  and  mission  time  since  the  differences  in 
values  would  be  well  within  the  intrinsic  measurement  error  levels  of 


made  to 
density  at 
correct  ion 
the  sensors. 


Step  4:  Calculation  of  che  Covariance  Matrix 

As  shown  in  Appendix  D,  the  covariance  matrix  of  the  difference  in  observa¬ 
tions  from  the  two  sensors  is  the  sum  of  the  measurement  error  covariance  matrices 
for  the  sensors.  Since  the  elements  of  the  measurement  error  covariance  matrix 
are  functions  of  range  and  azimuth,  they  should  be  recalculated  anytime  range  and 
azimuth  are  modified  in  the  bias  estimation  procedure.  There  are,  however,  two 
reasons  why  this  does  not  appear  to  be  worth  the  additional  burden  in  computer 
time.  The  first  is  based  on  the  assumption  that  the  biases  in  the  system  can  be 
kept  small  or,  at  least,  can  be  approximately  estimated  initially.  In  that  case, 
the  relative  change  in  matrix  elements  due  to  small  increments  in  range  and  azi¬ 
muth  will  be  correspondingly  small.  A  second  reason  stems  from  the  fact  that  the 
sensor  measurement  errors  are  also  involved  in  the  covariance  matrix  and  their 
values  arc  never  known  exactly. 

To  summarize  -  the  superpoints  are  first  corrected  for  known  o  estimated 
biases  and  then  adjusted  for  atmospheric  refraction  before  elements  of  the  dif¬ 
ference  covariance  matrix  are  computed. 

Step  5:  Function  Minimization 

Solving  for  the  bias  parameter  values  which  best  fit  a  given  set  of  super¬ 
points  from  each  of  two  sensors  is  a  problem  taken  from  maximum  likelihood  theory. 
As  developed  more  fully  in  Appendix  D,  the  solution  amounts  to  finding  the  minimum 
of  a  function  S  which  is  the  sum  of  the  squares  of  the  residual  vectors  scaled  to 
take  into  account  the  orientation  of  the  combined  measurement  error  ellipse.  The 
coordinate  transformations  from  range  and  azimuth  into  common  geodetic  coordinates 
cause  the  quantity  S  to  be  a  non-linear  function  of  the  biases. 

Several  approaches  were  considered  for  the  minimization.  One  possibility  is 
to  formulate  the  problem  into  the  familiar  matrix  equation. 
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A  x 
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where  x  is  the  vector  of  bias  components,  b  is  the  vector  of  observed  differences 
in  target  position  and  A  is  the  matrix  resulting  from  linearization  of  the  problem 
by  a  Taylor  expansion.  Subroutines  are  readily  available  for  finding  a  least- 
square  solution  to  this  problem  once  the  matrix  equation  is  set-up.  This  approach 
was  ultimately  abandoned  as  more  bias  components  were  introduced  into  the  model 
and  the  problem  of  developing  analytical  expressions  for  the  elements  of  the  A 
matrix  became  more  involved. 

Another  standard  approach  not  involving  the  calculation  of  derivatives  was 
also  investigated.  In  this  algorithm,  referred  to  as  the  grid  search,  the  func¬ 
tion  S  is  minimized  in  each  bias  parameter  separately.  It  offers  the  advantage  of 
straightforward  computer  programming  and  did  yield  useful  initial  results.  How¬ 
ever,  because  it  converges  very  slowly,  particularly  when  the  bias  parameters  are 
not  completely  independent,  it  was  not  deemed  suitable. 

A  third  method  was  applied  which  has  desirable  characteristics  of  rapid 
convergence  without  relying  on  analytical  derivations  or  the  evaluation  of  deriva¬ 
tives.  This  algorithm  was  proposed  by  Powell  (see  refs.  5  and  6),  and  is  based  on 
quadratic  convergence  and  some  properties  of  conjugate  vectors. 

Assuming  that  M  bias  components  are  to  be  solved  for,  the  algorithm  begins 
with  the  initial  best  estimate  of  the  bias  vector  (the  elements  of  which  are  the 
bias  components)  and  a  system  defined  by  the  M  linearly  independent  coordinate 
directions.  It  proceeds  as  two  loop  structures  (one  nested  within  the  other);  M 
iterations  of  the  inner  loop,  each  minimizing  S  results  in  the  establishment  of  a 
new  conjugate  vector.  The  function  S  is  then  minimized  along  the  new  conjugate 
direction  and  the  outer  loop  is  repeated  with  the  latest  estimate  of  the  bias 
vector  replacing  the  initial  best  guess. 

Each  minimization  of  S  along  a  conjugate  direction  is  accomplished  by  a  two- 
pass  increment  and  search  process.  In  the  first  pass  the  bias  vector  is  incre¬ 
mented  until  a  first  approximation  to  the  minimum  is  found.  In  the  second  pass 
the  increment  step  size  is  reduced  by  a  factor  of  10  and  the  search  is  repeated 
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until  a  new  minimum  is  detected.  As  a  last  step,  3-point  parabolic  interpolation 
is  employed  to  solve  tor  the  final  minimum  and  the  associated  value  of  the  bias 
vector  at  the  minimum. 

For  true  quadratic  functions  the  outer  loop  is  exercised  for  M  iterations 
before  the  minimum  is  reached.  This  constitutes  one  pass  through  the  algorithm 
and  is  used  as  the  basic  performance  unit. 

The  final  minimum  in  S  has  utility  as  a  measure  of  the  goodness  cf  fit  of  the 
bias  model  for  a  given  set  of  data.  S  is  a  statistical  quantity  having  a  chi- 
squared  probability  distribution  with  v  =  2N  -  M  -  1  degrees  of  freedom.  (N  is 
the  number  of  superpoints.)  Thus,  if  the  measurement  errors  in  range  and  asim>.  i.i- 
were  known  and  all  the  biases  accounted  for,  the  quantity  =  S/v  would  have  an 
expected  value  of  1.0  and  a  variance  of  2/v  (at  least  for  large  N,  see  Appendix  D, 
Section  E) . 


Finally,  even  though  S  is  not  a  true  quadratic  function,  the  Powell  algorithm 
was  found  to  converge  fairly  rapidly  toward  a  minimum.  This  is  demonstrated  in 
Figure  1  where  has  been  plotted  versus  pass  number.  After  only  one  pass  the 
value  of  is  se^n  to  decrease  from  a  value  of  3S.8  to  approximately  2.48. 


Step  6:  Estimating  the  Variances  of  the  Biases 


As  pointed  out  in  Appendix  D,  an  estimate  of  the  variance  of  the  bias  com¬ 
ponent  can  be  found  from  the  inverse  of  the  Hessian  or  curvature  matrix  H.  If 
is  taken  as  an  element  on  the  diagonal  of  the  inverse  of  H,  the  relation 


has  been  taken  as  a  first-order  approximation  to  the  variance  in  the  bias  estimate. 
Note  that  under  ideal  conditions  the  expected  value  of  S/v  is  approximately  1  and 
the  variance  would  depend  only  on  the  extent  of  cross  coupling  between  biases. 

VI.  EXPERIMENTAL  VERIFICATION  OF  ALGORITHM 

Data  recorded  simultaneously  from  two  sensors  was  used  to  exercise  the  self- 
registration  algorithm.  In  April  1977,  beacon  data  were  recorded  using  Air  Traffic 
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Control  Radar  Beacon  System  (ATCRBS)  interrogators.  One  sensor  was  located  at 
Lincoln  Laboratory,  Lexington,  Mass.,  and  the  other  at  1.1.  C.reen  Airport,  Provi¬ 
dence,  R.I.  Figures  2  and  3  are  photographs  of  the  two  sensors.  Each  ATCRBS 
interrogator  employed  a  monopulse  antenna  for  azimuth  angle  detcrminat ion.  Only 
four  or  five  interrogations  were  employed  as  the  antenna  scanned  by  the  target. 
Azimuth  was  determined  on  each  interrogation  and  the  results  were  averaged  to 
produce  a  very  accurate  azimuth  estimate.  The  sampling  rate  was  once  ever)'  4  to  5 
sec.  for  each  sensor.  Interrogations  alternated  between  Modes  A  (identity)  and  C 
(altitude).  The  identity  code  was  used  to  sort  out  The  replies  Trom  a  particular 
aircraft.  The  aircraft's  altitude  was  reported  from  a  barometric  altimeter  aboard 
each  aircraft  to  a  precision  of  100  ft.  Lata  was  recorded  simultaneously  from 
both  sensors  over  a  35-minute  period. 

As  described  in  Appendix  C,  the  separate  aircraft  reports  were  time  tagged. 
The  time  had  been  recorded  periodically  along  with  the  present  azimuth.  Target 
report  time  was  determined  by  interpolating  between  azimuth-time  reports.  Air¬ 
craft  altitude  reports  were  corrected  for  the  barometric  readings  at  sea  level. 
This  involved  adding  300  ft  to  all  altitude  reports. 

Figure  4  shows  the  paths  followed  by  the  aircraft  used  most  often  in  the 
analysis  which  follows.  The  marked  points  are  the  "superpoints".  Each  aircraft 
is  identified  by  its  code  number.  The  aircraft  tracks  are  plotted  in  a  Mercator 
proj  ection . 

After  creating  software  for  the  algorithm  as  outlined  in  the  last  section  the 
analysis  plan  proceeded  along  the  following  steps: 

1.  Verification  of  the  correctness  of  the  procedure  through 
the  use  of  simulated  target  reports. 

2.  Comparison  of  results  obtained  from  different  aircraft. 

3.  Investigation  of  some  approaches  toward  complete  automation 
of  the  sensor  sell -registration  problem. 

A.  Simulation 

It  became  apparent  after  processing  some  typical  target  data  that  the 
software  could  not  be  adequately  tested  using  real  measurements  alone.  Although 
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two  large  biases  were  known  to  exist  in  the  system,  little  could  be  surmised  about 
the  remaining  components. 

Simulated  data,  though  not  ideal,  does  at  least  provide  an  element  of  control 
and,  hence,  the  desired  ground  truth  for  general  testing  purposes.  In  the  in¬ 
terest  of  making  the  testing  realistic,  simulated  tracks  of  target  reports  were 
generated  from  real,  Lexington,  Mass.  data. 

Shown  in  Figure  S  are  the  steps  taken  in  going  from  actual  Lexington  target 
reports  to  a  file  of  simulated  superpoints  for  each  sensor.  The  RMS  measurement 
error  used  for  range  and  azimuth  was  that  assumed  for  the  real  data:  10  meters  and 
0.5  milliradian,  respectively.  The  systematic  (bias)  error  used  in  the  test  runs 
was  taken  as  representative  of  what  could  be  expected  in  ar.  operational  system. 
Notice  that  the  Providence  reports  would  agree  exactly  with  the  Lexington  reports 
if  the  noise  and  bias  additions  are  set  to  zero. 

Using  this  simulated  data,  two  different  algorithm  tests  were  performed.  The 
first  was  directed  at  the  question,  "Given  prior  knowledge  of  the  presence  of  a 
subset  of  bias  components  in  the  observations,  how  well  does  the  algorithm  do  in 
estimating  the  values  of  these  components0" 

To  answer  this  question,  Lexington  data  from  target  4655  was  used  to  generate 
simulated  data.  Five  biases  were  assumed  as  in  Table  II.  The  self-registrat ion 
algorithm  was  ap,  lied  with  the  results  shown  in  Figures  6  and  7  and  Table  II. 

TABLE  IJ 

SIMULATED  TARGET  REPORTS 


Known 

Computed 

A  2 1 1  (radj 

0.003 

0.00293 

R2  (n) 

150 

154 

A2I2  (rad) 

-0.002 

-0.0019S 

1 B1 2  (ms) 

-180 

-184 

ACT!)  (m) 

-100 

-no 

Assumed  Measurement  Error 

Range  RMS  -  10m 

Azimuth  RMS  -  0.0005  rad. 
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Fig.  5.  Generation  of  simulated  target  report  files. 
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Fig.  6.  Superpoints  (asterisks)  generated  from  simulated  Lexington 
data  and  error  vectors  between  sensors  magnified  200  times. 


Typical  plotted  output  from  one  of  these  test  runs  is  shown  in  Figure  6. 

Line  segments  in  the  diagram  are  "error"  vectors  drawn  from  Lexington's  observed 
target  positions  marked  by  asterisks  to  those  simulated  for  Providence.  The 
plot  is  a  Mercator  projection  and  for  illustration  the  error  vectors  are  scaled 
to  200  times  actual  length. 

After  only  one  pass  through  the  algorithm  the  results  in  Figure  7  were  obtained. 
The  error  vectors,  now  appreciably  reduced,  are  residuals  commensurate  with  the 
random  measurement  error  introduced  into  the  simulation. 

Displayed  in  Table  II  are  the  "known"  bias  components  as  compared  with 
values  computed  by  the  algorithm.  It  seems  reasonable  to  conclude  from  this  test 
that  given  a  priori  knowledge  of  the  existence  of  certain  bias  components,  the 
algorithm  is  capable  of  estimating  their  values  within  a  relative  error  of  a 
few  percent. 

The  second  level  of  algorithm  testing  addressed  the  more  realistic  situation 
in  which  the  presence  or  absence  of  individual  bias  components  is  not  known 
a  priori.  In  these  tests  all  biases  assumed  for  the  model  are  solved  for. 

Table  III  summarizes  the  output  obtained  from  one  of  these  runs  using  the 
same  simulated  input  data  of  Figure  6.  Also  tabulated  is  the  standard  deviation 
of  the  bias  estimate  computed  according  to  the  relationship  developed  in  Appendix 
D,  Section  C. 

Of  particular  note  in  the  table  are  bias  components  such  as  AZ21,  AZ31,  AZ22, 
AZ32  and  AZ42  which  have  "known"  values  of  zero  but  have  computed  values  signifi¬ 
cantly  different  from  zero.  As  a  consequence  of  solving  for  all  biases  in  the 
model,  the  "known"  components  with  non-zero  values  are  no  longer  as  accurately 
determined.  Some  of  the  discrepancy  in  reported  target  position  is  assigne 
the  other  components.  This  conclusion  is  also  reinforced  by  the  rather  high 
value  of  bias  standard  deviation.  For  example,  the  azimuth  offset  components 
AZ11  for  the  Lexington  sensor  and  AZ12  for  the  Providence  sensor  have  standard 
deviations  nearly  an  order  of  magnitude  greater  than  the  estimates  themselves. 
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TABLE  III 

BIAS  ESTIMATES  FOR  SIMULATED  DATA 

Aircraft  Code  =  46SS 

28  Superpoints 

Final  Value  of  =  .354 

COMPONENT  BIAS  ESTIMATE 

STANDARD  DEVIATION 

Lexington,  Mass. 

RC1 

-0.00083 

0.00698 

A21 1  (rad) 

0.00217 

0.033 

AZ21  (rad) 

-0.00154 

0.0024 

AZ31  (rad) 

-0.00141 

0.0026 

AZ41  (rad) 

-0.00064 

0.0026 

Providence,  R.I. 

RC2 

-0.00093 

0.0066 

R2  (m) 

167 

62 

AZ12  (rad) 

-0.00218 

0.032 

AZ22  (rad) 

-0.01366 

0.008 

AZ32  (rad) 

-0.00179 

0.0027 

AZ42  (rad) 

0.01141 

0.006 

TB1 2  (ms) 

-137.3 

55 

TB22 

0.00004 

0.00003 

LAT2  (rad) 

0.00001 

0.00011 

LON2  (rad) 

0.00002 

0.00056 

ACTD  (m) 

-112.6 

c;q 

For  comparison,  similar  results  using  real  data  are  presented  in  Table  IV. 
Despite  the  impressive  reduction  in  S  from  3S.8  down  to  2.4  which  is  a  measure  of 
the  goodness-ot -f it ,  the  high  values  of  standard  deviation  render  the  results 
somewhat  questionable.  In  fact,  when  other  aircraft  tracks  were  used,  quite 
different  values  for  the  biases  were  obtained. 


Thus,  further  insight  into  the  relationships  between  the  bias  components  was 
required  to  avoid  generating  highly  inaccurate  bias  estimates. 

To  generate  further  insight,  we  examined  the  Hessian  matrix  mentioned  breifly 
at  the  end  of  Section  V  and  introduced  during  the  development  of  the  bias  variance 
estimates  in  Appendix  D.  Because  of  its  relationship  to  the  curvature  of  S  in  the 
space  defined  by  the  biases,  the  Hessian  matrix  evaluated  near  the  minimum  in  $ 
can  reveal  the  degree  of  coupling  or  correlation  among  the  biases. 

Specifically,  if  h^  are  the  elements  of  the  Hessian  matrix  computed  by 
numerical  approximation  at  the  minimum  in  S,  then 


h. 


(h. 


h .  .  1 
3  )' 


1/2 


defines  a  coupling  coefficient  which  varies  between  -1  and  +1.  A  value  o£ 
a. .  near  1  means  that  an  increase  or  decrease  of  the  l  bias  component  has 


ij 
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the  same  effect  on  S  as  a  corresponding  change  in  the  j  bias  component. 


Table  V  was  prepared  for  the  aircraft  having  an  ID  code  of  4655.  It  is 
quite  apparent  that  several  biases  are  tightly  coupled  for  this  set  of  superpoints. 
Notably  high  coefficients  are  found  for  the  pairs:  (R2,  RC2) ,  (LAT2,  ACT'D) , 

(A212,  AZ22),  (AZ21,  AZ4L)  and  AZ22,  AZ42).  The  results  using  other  sets  of 
superpoints  are,  in  general,  different  although  for  certain  bias  pairs  such  as 
(R2 ,  RC2) ,  (A21 1 ,  AZ21) ,  (LAT2,  ACTD)  and  (TB12,  TB22)  the  are  consistently 
higli  and  positive. 

Thus,  depending  on  the  spatial  distribution  of  an  aircraft's  superpoints,  it 
may  not  always  be  possible  to  unambiguously  estimate  those  biases  contributing  to 
systematic  error  in  the  data.  Better  results  will  be  obtained  by  not  solving  for 
both  members  of  a  highly  coupled  pair.  A  similar  result  has  been  reported  by 
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Lexington,  Mass. 


Providence,  R.I, 


TABLE  IV 

BjiAS  ESTIMATES  FOR  REAL  DATA 


Aircraft  Code  =  46SS 
23  Superpoints 

Assumed  Bias:  AZ12  =  -.2410  (rad) 


Final 

Value  of  S  =24 

v 

component 

BIAS  ESTIMATE 

STANDARD  DEV I A 

RC1 

0.00469 

0.0144 

AZil  (rad) 

-0.  00344 

0.0320 

AZ2 1  (rad) 

-0.00163 

0.00730 

AZ31  (rad) 

-0. 00246 

0.00701 

AZ41  (rad 

-0.0026S 

0.00767 

RC2 

-0. 00568 

0.0132 

R2  (m) 

-425 

195 

AZ12  (rad) 

-.244 

0.0308 

AZ22  (rad) 

0.008 

0.0217 

AZ32  (rad) 

-0.0245 

0.00984 

AZ42  (rad) 

-0.00331 

0.0180 

TB1 2  (ms) 

-124 

174 

TB22 

-0.0000447 

0.000111 

LAT2  (rad) 

-0.0000344 

0.000193 

L0N2  (rad) 

-0.0000377 

0.000543 

ACTD  (m) 

76.7 

141 
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TABLE  V 

COUPLING  COEFFICIENTS  >0.5 


Aircraft  Code  =  4655 


A21 1 

LAT2 

LON  2 

.69 

AZ3 1 

.60 

AZ3 1 

.83 

ACTD 

.96 

AZ4 1 

-.70 

RC1 

LON  2 

LAT2 

.49 

AZ21 

.69 

AZ32 

-.73 

AZ41 

-.65 

ACTD 

.52 

AZ22 

.73 

AZ42 

.69 

R2 

TBI  2 

RC2 

.96 

TB22 

.93 

LAT2 

.86 

AZ3 1 

.69 

ACTD 

.83 

AZ1 2 

AZ21 

LON2 

.72 

AZ41 

-.94 

AZ22 

.93 

AZ42 

.84 

RC2 

AZ22 

LAT2 

.76 

AZ42 

.97 

AZ31 

.75 

ACTD 

.69 
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Oliver  (ref.  7)  when  using  a  least -squares  procedure  to  determine  the  biases 
in  the  angular  pointing  of  an  antenna. 

You  will  recall  that  in  Table  II  only  a  few  biases  were  solved  for  and 
that  the  predicted  biases  were  quite  close  to  the  known  values,  f-'or  those  super¬ 
points  the  highest  coupling  coefficient  was  0.79  for  the  pair  IR2,  ACTO) . 

B.  Applying  Bias  Estimates  to  Several  Aircraft 

In  order  to  generate  bias  estimates  for  application  to  other  aircraft  the 
superpoints  shown  in  Figure  4  were  prepared.  Three  aircraft  having  identification 
codes  of  4655,  4543  and  3430  were  individually  processed  and  then  combined  into 
one  large  data  set  of  64  superpoints.  As  can  be  seen  from  Figure  4,  these 
superpoints  have  the  reasonably  broad  geographical  distribution  desired  for  bias 
estimation . 

Next,  a  set  of  eight  bias  components  was  somewhat  arbitrarily  selected. 

These  were  to  be  the  only  components  solved  for  and  it  was  required  that  they  be 
only  weakly  coupled.  The  selected  components  and  corresponding  coupling  coef¬ 
ficients  are  given  in  Tame  VI.  Some  items  require  special  note:  (1)  the  com¬ 
ponents  ACT1 ,  ACT2  and  ACT3  refer  to  the  transponder  delay  biases  for  aircraft 
with  codes  of  4655,  4543  and  3430,  respectively,  (2)  the  matrix  is  symmetric  by 
definition  and  only  the  upper  triangular  portion  is  shown,  and  (3)  the  biases  are 
only  weakly  coupled  as  desired  with  the  highest  coefficient  being  only  0.69  for 
the  pair  LAT2 ,  ACT2. 

The  bias  estimates  and  standard  deviations  calculated  for  the  64-superpoint 
set  are  listed  in  Table  VII.  As  expected,  and  due  primarily  to  the  restriction  on 
the  degree  of  coupling  permitted,  the  bias  standard  deviations  are  s len  to  be 
acceptably  small. 

The  resulting  value  of  (3.10)  gives  an  estimate  of  the  variances  of  any 
residual  differences  in  the  data  points  after  application  of  the  estimated  biases. 
If  =  1.76  is  multiplied  by  the  assumed  range  and  azimuth  errors  (10  m  and  0.5 
mrad) ,  the  results  17  m  and  0.88  mrad  represent  the  residual  differences  between 
data  points. 
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TABLE  Vil 

BIAS  ESTIMATES  FROM  64-SUPERPOINT  SET 


BIAS  ESTIMATE 

STANDARD  DEVIATION 

AZ11  (rad) 

-0.000693 

0.000126 

R2  (m) 

-195 

9.1 

AZ1 2  (rad) 

-.2409 

0 . 0001  S3 

LAT2  (rad) 

-0.0000191 

0.0000019 

TB12  (msec) 

-194.5 

34.1 

ACT1  (m) 

1SS.  s 

10.8 

ACT2  (in) 

123. 1 

10.5 

ACT3  (m) 

249.3 

11  .4 

ASSUMED  INITIAL  BIAS:  AZ12  *  -.2410 
INITIAL  Sv  =  28.1  Sv  (AFTER  TWO  ITERATIONS)  =  3.10 
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Having  established  a  set  of  bias  estimates,  the  next  point  of  interest  was  to 
investigate  how  well  systematic  error  in  target  reports  from  other  aircraft  could 
be  accounted  for. 

For  the  following  reasons  it  was  decided  to  omit  the  time  offset  lias  T812. 

In  order  to  correctly  process  report  times  from  the  Providence  sensor,  the 
additional  preprocessing  described  in  Appendix  C  would  have  been  required  for 
each  aircraft.  Moreover,  even  if  the  preprocessing  were  performed  and  the  time 
bias  applied,  the  resulting  adjustments  would  only  be  on  the  order  of  40-50  meters 
and  always  in  the  direction  of  aircraft  motion.  Thus,  at  the  plotting  scales  em¬ 
ployed,  the  combined  tracks  would  appear  about  the  same  whether  or  not  a  time  bias 
was  applied  to  the  Providence  reports. 

Transponder  delay  is  a  bias  which  is  different  for  each  aircraft  and  not 
known  b  forehand.  The  approach  taken  here  in  correcting  the  reports  from  other 
aircraft  was  to  use  the  mean  value  of  the  three  transponder  delays  which  had  been 
calculated  and  are  shown  in  Table  VI  i. 

When  the  corrected  target  reports  are  plotted  individually  for  each  sensor, 
one  finds  the  maximum  separation  of  tracks  occurs  when  the  aircraft  is  flying  in 
the  airspace  between  the  two  sensors  and  at  a  heading  cf  90°  with  respect  to  the 
line  connecting  them.  In  fact,  rough  estimates  of  the  transponder  delay  bias  can 
be  estimated  using  this  data. 

To  recount  the  discussion  above,  target  reports  for  additional  aircraft  were 
corrected  tor  the  biases  designated  AZ11,  R2,  AZ12  and  LAT2  as  well  as  a  number 
representing  an  average  transponder  delay  error.  An  independent  set  of  throe 
aircraft  were  used  to  estimate  the  bias  values  which  are  listed  in  Table  VII. 

Although  seven  additional  aircraft  were  considered,  only  a  few  typical  cases 
will  be  presented.  As  was  true  earlier,  a  Mci  cator  projection  is  employed. 

Shown  in  Figure  8  for  aircraft  code  3502  are  the  tracks  as  they  are  reported 
at  Lexington  and  Providence  before  any  biases  are  accounted  for.  The  large 
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azimuth  offset  of  approximately  14°  in  the  Providence  measurements  is  quite  con¬ 
spicuous.  After  bias  removal  the  reports  form  the  common  track  of  Figure  9.  The 
value  of  dual  sensor  coverage  is  demonstrated  in  Figure  10  where  the  gap  and  low 
quality  of  the  Lexington  reports  is  filled  in  quite  adequately  by  the  supporting 
Providence  data. 

Figures  11  and  12  are  illustrations  of  the  aircraft  with  code  2417,  Cne  in¬ 
teresting  feature  here  is  the  apparent  track  segment  at  an  azimuth  of  around 
130°  from  Lexington  which  is  reported  only  by  Lexington  and  not  Providence.  This 
segment  was  found  to  be  a  false  track,  i.e.,  a  track  due  to  reflections  from  an 
obstruction  (a  large  semicircular  antenna  reflector)  in  the  immediate  vicinity  of 
the  Lexington  site. 

This  porticn  of  the  coverage,  expanded  in  Figure  13,  also  shows  how  the  Lex¬ 
ington  reports  can  be  used  in  lieu  of  those  from  Providence  just  before  the  sharp 
bend  in  the  track.  The  nearly  one-kilometer  scatter  in  the  positions  of  the 
Providence  reports  is  apparently  due  to  low-eievation-angle  propagation  anomalies. 

C.  Automatic  Bias  £st irnatiop. 

In  this  section  two  alternatives  for  completely  automating  Lhe  bias  estima¬ 
tion  are  discussed.  It  was  shown  in  the  last  section  that  fairly  good  results  can 
be  obtained  by  accounting  foi  only  four  bias  components.  What  is  desired  now  is  a 
procedure  which  both  chooses  the  bias  components  given  a  set  of  superpoints  and 
then  solves  for  them.  The  objective  is  to  find  that  combination  of  biases  which 
yields  a  low  value  of  Sy  and  at  the  same  time  has  consistently  small  coupling 
coefficients . 

We  consider  two  procedures.  Procedure  A  outlined  below  places  primary  em¬ 
phasis  on  the  value  of  S  and  secondary  emphasis  on  the  degree  of  coupling  between 
the  components.  In  Procedure  B,  just  the  reverse  is  true.  In  both  it  is  assumed 
that  the  biases  in  one  of  the  two  sensors  (the  Lexington  sensor)  are  small  enough 
so  that  the  elements  of  the  Hessian  matrix  can  be  estimated  to  within  a  reasonable 
degree  of  accuiacy.  This  is  accomplished  using  simulated  reports  from  the  other 
sensor  exactly  in  the  manner  described  in  the  earlier  section  on  simulation,  but 
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Aircraft  code  3502  after  bias  removal 
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Fig.  10.  Expanded  view  of  portion  of  aircraft  code  3502  tracks 
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Fig.  H,  Aircraft  code  2417  as  seen  by  two  sensors. 
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Fig.  13.  Expanded  view  of  aircraft  code  2417.  Observe  spurious 
reflections  from  large  obstruction  at  about  130°  azimuth. 


without  the  addition  of  biases.  Comparisons  with  results  obtained  from  real 
measurements  have  shown  that  this  method  works  quite  well. 

Procedure  A  begins  by  computing  Sv  when  each  of  K  biases  is  solved  for  in¬ 
dividually.  Recall  that  for  one  aircraft  K  is  equal  to  16. 

s'-  ’■ 

The  first  bias  selected  is  that  which  results  in  the  lowest  value  of  S  . 

v 

This  process  is  then  repeated  K-l  times  with  two  biases  (the  one  already  selected 
on  the  previous  pass  and  one  candidate  from  the  list  of  remaining  components). 

The  candidate  bias,  i.e.,  the  one  for  which  is  a  minimum,  is  added  to  the 
selected  list  only  if  the  standard  deviations  of  the  already  selected  biases  are 
not  greater  than  some  prescribed  multiple  of  their  standard  deviation  when  first 
selected. 

Thus,  biases  which  might  have  a  negative  effect  on  those  already  selected 
components  because  of  tight  coupling  are  discarded  and  not  considered  further. 
Processing  terminates  when  each  of  the  K  bias  components  has  been  either  selected 
or  discarded. 

Table  VIII  summarizes  a  run  of  procedure  A  using  23  superpoints  and  a  mul¬ 
tiple  of  1.5  as  the  selection  parameter.  The  order  of  entries  in  the  table  re¬ 
flects  the  order  in  which  the  biases  were  selected.  Thus,  AZ12  was  the  largest 
bias  found  even  when  an  initial  value  of  -.2300  radian  (m.- 13.2  deg)  was  assumed 
for  it. 

One  variation  on  this  method  would  be  to  use  the  degree  of  coupling  between 
the  candidate  bias  and  the  previously  selected  biases  as  an  accept/reject  cri¬ 
terion.  This  would  reduce  the  amount  of  computation  since  the  full  Hessian  matrix 
is  computed  only  once  at  the  beginning  of  the  run  and  the  matrix  inversion  per¬ 
formed  in  computing  the  bias  standard  deviation  could  be  avoided. 

Procedure  B  places  greater  emphasis  on  the  degree  of  coupling  between  the 
bias  components  and  then  utilizes  the  calculated  S,  values  as  a  final  selector. 

As  was  the  case  in  Procedure  A,  the  elements  of  the  Hessian  matrix  are  cal¬ 
culated  and  from  them  the  coupling  coefficients  between  the  K  bias  components. 
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COMPONENT 

AZ12 

AZ31 

RC1 

TB27 

LON  2 

R2 

AZ32 

AZ41 

AZ42 


TABLE  VIII 

AUTOMATIC  BIAS  SELECTION  AND  ESTIMATION 
PROCEDURE  A 
23  Superpoints 
Aircraft  Cede  =  46SS 
Assumed  Biases:  AZ12  ■  -.2300 
Multiplier  =1.5 
Final  Value  of  =»  4.01 


ESTIMATED  VALUE 

-  .2405 
-0.0089 
0.0019 
-0.000138 
-0.000013 
-35.74 
0.00460 
-0.0014 
-0.00081 


STANDARD  DEVIATION 

0.00053 

0.00267 

0.000301 

0.000045 

0.000006 

17.99 

0.0046 

0.0030 

0.0039 
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Then,  establishing  a  cut-off  on  the  magnitude  of  the  coupling  permitted,  N  com¬ 
binations  of  J  biases  each  are  formed.  These  N  combinations  each  constitute  a 
weakly  coupled  set  such  that  no  set  with  greater  than  J  biases  could  be  created. 

For  each  of  the  N  combinations  is  calculated,  and  the  final  choice  is 

based  on  which  results  in  minimum  S^.  In  cases  where  several  combinations  yield 
values  of  S^,  all  within  the  calculated  standard  deviation  of  S^,  approximately 
2/v,  then  the  combination  with  the  uniformly  lowest  coupling  would  be  selected. 

The  64-superpoint  data  set  described  earlier  and  shown  in  Figure  4  was  the 
basis  for  testing  Procedure  B.  When  a  limit  of  0.5  was  placed  on  the  largest 
coupling  coefficient  it  was  found  that  S4  combinations  of  eight  biases  each  could 
be  formed.  There  were  none  with  nine  or  more  combinations.  Table  IX  lists  values 

calculated  for  the  combination  which  resulted  in  minimum  S  . 

v 

Both  procedures  have  attributes  which  are  not  satisfactory. 

1.  They  require  selection  parameters  which  can  affect  the  final 
outcome. 

2.  They  are  dependent  on  initial,  assumed  bias  estimates. 

3.  They  consume  a  great  deal  of  computer  processing  time. 

Of  the  two  approaches  Procedure  B  seems  to  be  preferable.  It  is  faster  and 
relies  less  on  assumptions  concerning  initial  bias  values.  Moreover,  intuition 
seems  to  favor  the  generation  of  conservative  estimates  of  only  a  few  biases  at 
the  expense  of  having  to  accept  larger  residual  errors.  Then  as  more  data  becomes 
available  in  additional  coverage  areas,  new  bias  estimates  can  be  optimally  com¬ 
bined  with  previous  estimates  in  the  manner  suggested  in  Appendix  D. 

VH.  CONCLUSIONS 

The  problem  of  accurately  registering  sensors  in  a  network  has  been  studied. 
There  are  many  possible  sources  of  calibration  errors  affecting  the  registration 
accuracy.  These  include  such  things  as  azimuth,  range  and  time  offset  biases, 
errors  in  sensor  location,  etc.  When  all  of  these  error  sources  are  eliminated, 
and  when  accurate  coordinate  conversion  equations  are  employed,  various  sensor 
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TABLE  IX 

AUTOMATIC  BIAS  SELECTION  AND  ESTIMATION 
PROCEDURE  B 
64  Supcrpoints 

Assumed  Bias:  AZi2  =  -.2400 
Coupling  Coefficient  Limit  =  0.5 
Sv-  3.16 


COMPONENT 

ESTIMATED  VALUE 

STANDARD  DEVIATION 

AZ1 1  (rad) 

-O.OOOS4 

0.00013 

RC1 

0.00019 

0.000057 

RC2  (m) 

-88.7 

6.2 

AZ31  (rad) 

-0.0129 

0.00135 

AZ22  (rad) 

-0.0049 

0.00085 

TB22 

-0.00011 

0.0000256 

ACT1  (m) 

57. 7 

6.15 

ACT3  (m) 

148.1 

6.97 
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outputs  will  register  with  the  accuracies  determined  by  measurement  accuracies  cf 
the  individual  sensors. 

A  self-registration  algorithm  was  developed  to  check  for  bias  errors  and  to 
find  a  set  of  bias  corrections  which  will  cause  good  registration  over  the  whole 
common  coverage  of  two  sensors.  Usually  ail  of  the  possible  biases  cannot  be 
determined  unambiguously  due  to  the  coupling  effect  between  biases.  This  coupling 
effect  and  its  impact  on  the  variance  of  the  bias  estimates  are  also  evaluated  by 
the  self-registration  algorithm  and  can  be  used  as  criteria  for  the  choice  of  a 
subset  of  all  possible  biases  so  as  to  provide  adequate  registration  accuracy. 

The  self-registration  algorithm  was  tested  by  analyzing  both  simulated  and 
real  data  from  two  sensors  spaced  about  45  nmi  apart.  The  results  demonstrated 
the  essential  accuracy  of  the  algorithm  and  provided  final  rms  registration 
accuracies  of  about  one  milliradian  in  azimuth  and  20  meters  in  range  which  is 
about  twice  the  measurement  accuracies  of  the  individual  sensors.  Indeed  the 
registration  was  good  enough  to  easily  discern  areas  of  increased  sensor  error. 

This  increased  error  was  generally  confined  to  low-elevation  aircraft  and  could 
usually  be  traced  to  obstructions  near  the  sensor. 

The  self-registration  algorithm  developed  in  this  study  is  not  considered  a 
cure-all  for  registration  problems.  The  factors  limiting  registration  accuracy 
should  be  corrected  at  their  source.  Tilt  meters,  north-seeking  gyros  and  ex¬ 
ternal  sensor  posi t ion -location  systems  should  be  employed.  Accurate  clocks  and  a 
system  for  clock  registration  should  be  used.  Range  offsets  and  antenna  tilt 
should  be  carefully  calibrated.  Because  of  the  coupling  among  the  biases  de¬ 
scribed  above,  the  self-registration  algorithm  should  only  be  used  as  an  overall 
check  of  registration  accuracy  to  discover  items  not  properly  calibrated  or  to 
refine  a  few  biases. 

The  ability  to  accurately  register  the  output  of  sensors  with  overlapping 
coverage  now  allows  the  development  of  multi-sensor  trackers.  A  multi -sensor 
tracker  accepts  target  reports  from  more  than  one  sensor  and  forms  a  single  track 
for  each  aircraft.  Multi-sensor  tracking  provides  the  advantages  enumerated  in 
the  Introduction;  namely,  improved  track  quality,  identity  maintenance,  large  area 
coverage  including  low  flyers  and  resistance  to  enemy  countermeasures.  In  addition, 
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a  mulct  censor  network  provides  the  basis  for  inclusion  of  other  typos  of  in¬ 
formation  such  as  radio  direction  finder  reports  and  aircraft  position  reports 
derived  from  on-board  aircraft  navigation  systems. 

The  registration  accuracy  demonstrated  in  this  study  is  considerably  better 
than  that  customarily  achieved.  In  order  to  take  full  advantage  of  this  increased 
accuracy  the  transmission  of  more  significant  bits  in  the  reporting  message  for¬ 
mats  will  be  required.  The  formats  should  also  include  target  report  time  tags. 
Lastly,  it  is  recommended  that  all  reports  be  in  geodetic  coordinates  to  allow 
easy  correlation  with  other  data  already  in  geodetic  coordinates  and  to  allow  an 
integrated  system  of  sensors  to  grow  naturally  over  an  extended  area. 
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APPENDIX  A 


PROPAGATION  DELAY  IN  ATMOSPHERE 


The  correction  to  range  for  the  effect  of  atmospheric  refraction  is  based 
on  the  following  steps: 


1.  Calculate  target  elevation  angle  n  from  the  expression: 


.  -1 

n  =  sin 


H2  -  H2  *  2  E  (H  -  H  ) 
r  v  r 

2  R  (E  +  H  ) 
r 


R 


2 


0 


where: 


H  =  target  altitude 

Hr  =  radar  altitude 

R  =  target  range 

E  =  mean  earth  radius 

2.  Kith  range  and  elevation,  refer  to  Table  A-l  for  the  correction  value. 

a.  Use  logarithmic  interpolation  for  ranges  and  elevations  not 
found  directly  as  entries  into  the  Table. 
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APPENDIX  B 


COORDINATE  TRANSFORMATIONS 


A.  Cartesian  Coordinates  Given  Geodetic  Coordinates 

Ke  first  derive  equations  for  converting  between  geodetic  coordinates  lati¬ 
tude  <J>,  longitude  and  altitude  li,  above  mean  sea  level  and  x,  v,  z  coordinates 
whose  center  is  at  the  earth's  center.  As  shown  in  Figure  B-l,  the  z  axis  co¬ 
incides  with  the  earth's  axis  of  rotation  and  the  x  axis  passes  through  the 
Greenwich  meridian.  The  eartii  has  a  radius  "a"  at  the  equator  and  "b"  at  the 
poles.  Latitude  is  defined  as  the  angle  p  that  a  normal  at  the  earth's  surface 
makes  with  its  projection  on  the  equatorial  plane.  The  eccentricity  "e"  is  defined 

as:  ,2.2 

z  a  -  o 

e  =  - ^ -  (Bl) 

a 


From  the  properties  of  an  ellipse  (see  Figure  B-2)  it  can  be  shown  that  the 
normal  nf.  terminating  at  the  major  axis  equals: 


a  (1 


<J>) 


1/2 


(B2) 


The  normal  nm  produced  to  the  minor  axis  equals: 


N  = 


(1  -  e2  sin2  p)1/2 


(B3) 


Using  these  values  and  referring  to  Figure  B2,  it  is  easily  seen  that  the  follow¬ 
ing  equations  convert  geodetic  coordinates  to  x,  y,  z  coordinates. 


a  (1  -  e  ) 


,,  2.2  .,1/2 

(1  -  e  sin  p) 


H 


sin  P 


(B4) 


Letting  line  segment  cd  be  called  s  in  Figure  B-2, 
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MERIDIAN 

Fig.  B-l .  Coordinate  geometries. 


Fig.  B-2.  Cross  section  of  Earth  through  its  axis 
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A 


M  2a  .  2  1/2  +  H 

,(1  -  e  sm  <£)  J 

cos  <J> 

(B5) 

s  cos  A 

( B6) 

s  sin  X 

(B7) 

B .  Radar  Coordinates  to  Geodetic  Coordinates 

Accurate  conversion  from  radar  coordinates  to  geodetic  coordinates  can  be 
accomplished  in  three  steps.  The  radar  coordinates  [range,  azimuth,  and  elevation) 
are  first  converted  into  local  Cartesian  coordinates.  The  local  coordinates  are 
then  rotated  and  translated  into  earth-centered  Cartesian  coordinates  finally, 
the  earth-centered  coordinates  are  converted  to  geodetic  coordinates  (latitude, 
longitude,  and  altitude  above  mean  sea  level).  We  will  describe  each  of  these 
steps  in  turn 

(1)  The  local  Cartesian  coordinates  are: 

x^  =  k  cos  r,  sin  6  (B8) 

=  R  cos  n  ccs  G  (B9) 

z^  =  R  sin  n  (BIO) 


where 

A 

points 

toward  the  zenith, 

points 

toward  the  north. 

X£ 

points 

toward  the  east, 

e 

is  the 

azimuth  measured  clockwise  from  north 

R 

is  the 

slant  range,  and 

n 

j  s  the 

elevation  angle  to  the  target. 

It  is  assumed  that  !\  and  n  have  been  corrected  for  propagation  effects  in  the 
earth's  atmosphere. 

(2)  The  local  coordinates  are  next  transformed  into  earth-centered  coordinates. 
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«uu 


(Bll) 


where  is  the  rotation  matrix  which  rotates  the  local  coordinate  system  to  align 
it  with  the  earth-centered  coordinate  system. 


-sin  X  -cos  A  sin  <6 

r  r  r 


cos  A  -sin  A  sin  6 

r  r  r 


cos 


cos  A  cos  d> 
r  r 


sin  X  cos  & 
r  r 


sin 


(B12) 


The  second  vector  in  the  above  equation  is  the  location  of  the  center  of  the  radar 
in  earth-centered  coordinates.  This  is  determined  using  Eqs.  (B4)  through  (37)  and 
using  <p,  A,  and  H  equal  to  the  coordinates  of  the  radar. 

(3)  Finally,  we  convert  the  earth-centered  Cartesian  coordinates  to  geodetic 
coordinates.  To  solve  for  the  latitude  <p  and  altitude  H,  it  would  be  logical  to 
invert  equations  (B4)  and  (B5) ,  but  this  would  necessitate  the  solution  of  a 
higher  order  equation.  Instead,  we  first  solve  for  H  using  the  original  data 
and  then  use  H  in  equations  (B4)  and  (B5)  to  solve  for  <p. 

Using  the  construction  of  Figure  B-3  and  the  law  of  cosines,  we  find 


7  7  l/? 

H  =  {(E  «•  Hrr  -►  R  +  2R  (E  ♦  Hr)  sin  n)  J 17  ^  -E 


(B13) 


where:  H  is  the  radar's  altitude  and  E  is  the  mean  earth  radius,  (a  <  b)/2. 
r 

Using  the  mean  earth  radius  instead  of  solving  exactly  in  spheroidal  coordinates 
introduces  a  maximum  error  of  25  m  in  H  at  a  range  of  300  km  increasing  as  R2. 
This  is  a  negligible  height  error  and,  in  addition,  causes  a  still  smaller  error 
in  determining  latitude  <$. 

Next  in  determining  <p,  the  square  root  term  is  eliminated  between  (B4)  and 
(B5),  yielding: 
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tan  $  = 


p  (1  -  e  )  +  e2  H  cos<j> 


(B14) 


,  ,  2  2.1/2 
where:  p  =  (x  +  y  )  . 

This  can  be  evaluated  ir.  two  steps.  A  first  estimate  of  4>  is  made: 


tan  = 


p  (1  -  e  ) 


(B15) 


Then  the  final  value  is  determined 


tan  $  = 


2  2 

p(l-e)+e  H  cos  0 


(B16) 


The  size  of  the  second  term  in  the  denominator  of  (B14)  is  about  1/3000  of 
the  first  term  for  a  high-altitude  aircraft.  Thus,  the  approximation  of  Eq.  (B16) 

7 

is  correct  to  about  one  part  in  10  ,  which  is  an  entirely  negligible  error. 
Finally,  the  longitude  is  determined: 


tan  X  =  £ 

x 


( B 1 7  ) 


The  above  conversion  process  might  seem  rather  lengthy.  It  actually  in¬ 
volves  the  determination  of  only  six  trigonometric  functions,  24  multiplies  or 
divides  and  two  square  roots. 

C .  Beacon  Coordinates  to  Geodetic  Coordinates 

A  beacon  transponder  system  produces  height  instead  of  elevation  angle.  Using 
Eq.  (B13J, 


z„  -  R  sin  n  - 


H2  -  H  2  +  2  (H  -  H  )  E  -  R2 
r  r 

2  (E  +  H  ) 
r ' 


(B18) 


R  1  -  (— )2  sin  6 


(B19) 


y£  =  R  h  -  (j/-)2  cos  e 


( B20) 
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equations  ( B 1 1 )  to  (B17)  are  then  used  to  find  the  earth  coordinates. 

D.  Geodetic  Coordinates  to  Local  Coordinates 

It  will  often  be  required  to  transform  from  geodetic  coordinates  (latitude, 
longitude,  and  altitude)  to  local  coordinates.  For  instance,  local  coordinates 
(azimuth,  elevation,  and  range)  are  required  to  aim  an  air  defense  missile  against 
a  target. 

(1)  Geodetic  coordinates  ($,  A,  H)  are  first  converted  to  earth-centered 
(z,  y,  z)  using  Eqs.  (B4)  through  (B7). 

(2)  Local  coordinates  (x^,  >^,  z are  generated  using: 


(B21) 


T 

where:  x_ ,  y  ,  2^  are  the  coordinates  ot  the  air  defense  weapon  and  T  is  the 
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~  T 
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X 

y£ 

=  T 
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-T  1 

R 

>'r 

_z_ 

■7 

u  r J 

transpose  of  the  rotation  matrix  in  Eq.  (B12). 
Finally,  R,  G  and  n  are  determined. 

R  =  yl*  : 

X£ 

tan  0  =  — 


2,1/2 
0  > 


(B22) 

(B23) 


tan  n  = 


[2 

\l*g  + 


(324) 


y£ 


E .  Reference  Ellipsoids 

The  earth  model  employed  in  this  investigation  is  the  DoD  World  Geodetic 
System  1972  reference  ellipsoid.  Another  model,  often  used  in  local  surveying, 
is  called  the  Clarke  1866  ellipsoid.  The  important  parameters  of  both  models 
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are  given  in  Table  B-l.  Bias  estimation  test  runs  indicated  little  sensitivity 
to  choice  of  model  at  the  latitudes  and  longitudes  considered  here. 

TABLE  B-l 

REFERENCE  ELLIPSOIDS 


SEMI-AXES 

MODEL  MAJOR  A _ MINOR  B _ FLATTENING  FACTOR 

DoD  WGS-72  637813S  m  63567S]  m  1/298.26 

Clarke  1866  6378206  m  63S6S8S  m  1/29S.0 


F.  Single  Precision  Coordinate  Transformation 

The  exact  coordinate  transformations  described  above  require  double 
precision  aritlunetic. 


A  very  accurate,  single  precision  coordinate  transformation  from  sensor 
coordinates  to  geodetic  coordinates  is  obtained  by  expansion  of  the  earth- 
surface  arc  length  in  series  form  (ref,  9).  In  what  follows  P, ,  1!,  and  0  axe 
the  target  range  and  height  above  sea  level  and  azimuth.  Hr>  <p  ,  and  X ^  are 
the  radar  height  above  sea  level  and  latitude  and  longitude,  and  a  and  e2  are 
the  earth's  radius  at  the  equator  and  eccentricity  squared.  For  sensor 
ranges  less  than  350  km  the  central  angle  represented  by  the  arc  length  c  is 
given  to  one-meter  accuracy  by: 

3 

n,  °o 

o  =  2  arc  sin  ( o  / 2)  -  a  *  ( B 2 5 ) 

o  o  24 


where: 


and 


(B26) 


(B27) 
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The  target's  latitude  is  given  by  first  calculating: 


tan  P1  2  2 

<5  d>  =  o  cos  0  -  I - t —  •  — |  o  sin  0 

2  p2 

(1  +  3  tan"  tp  o  2 
-  [ - - - -  (_)  ]  o  sin^  0  cos  O 

p  ^ 

1(1  +  3  tan^  <J  -  (9  +  IS  tan^  $.)  cos^  OJ  tan  <p  o4  sin"  6  (B28) 

r  l  yj  2  r 


-  3  e“  sin  4>  cos  <p 

<p  =  <*>  +  6$  -  (6$)"  ( - — — — 2 - ^  (b2S) 

2  (1  -  e‘  sin  (j)  ) 

r 

The  above  equation  for  target  latitude  was  derived  using  the  law  of  cosines 
for  spherical  triangles  which  gives  cos  o  and  expanding  in  a  series  about  the 
sensor  latitude.  Using  the  law  of  sines  for  spherical  triangles,  we  find: 


sin  0  sin  o. 

A  -  A  +  arc  sin  ( - 2 - -) 

r  cos  $  ; 

0.  sin  0  a?  sin  9  sin-5  0 

X  +  ^  ^  ^  i- 

r  cos  <6  *  6  cos  6  ~  T~, 

6  cos 


where: 


2  2 

o  P  j  (1  -  e  Sin  <p  ) 


The  above  equations  were  checked  for  accuracy  by  calculating  latitudes  and 
longitudes  at  many  azimuth  angles  for  sensor  ranges  out  to  350  km.  The  results 
were  compared  with  those  calculated  using  the  very  accurate  equations  presented 
in  the  first  part  of  this  Appendix.  The  maximum  differences  are  presented  in 
Figs.  B  4  and  8-S.  We  see  in  Fig.  B-4  that  adequate  accuracies  are  achieved  up 
to  about  45  degrees  latitude  using  a  four-term  expansion.  This  accuracy  can 
be  extended  to  70°  by  using  a  fifth  term  (the  fourth  term  in  Eq.  B28).  These 
curves  are  independent  of  target  height  to  well  above  18  km. 


61 


E  1,000 1 


TARGET  range  •  350km 
ALTITUDE  •  1 0  km 
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RADAR  LATITUDE  ld*g) 

Fig.  B-4.  Maximum  errors  introduced  by  using  series  approximations 
for  target  latitude. 
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APPENDIX  C 

DESCRIPTION  OF  PRE-PROCESSING 


A.  Introduction 

Data  employed  in  this  study  were  recorded  simultaneously  at  each  of  two 
Air  Traffic  Control  Radar  Beacon  System  (ATCRBS)  sites,  Lexington,  Mass,  and 
Providence,  R.I. 

The  target  reports  which  comprise  the  basic  data  have  already  been 
preprocessed  by  reply  correlation  and  centroiding  into  target  reports  as  well 
as  target-to-track  a~.  jciation.  An  extensive  documentation  of  that  surveillance 
processing  software  is  available  in  ref.  11  and  it  should  be  consulted  for 
additional  detail. 

In  this  study,  two  data  sources  have  been  used  to  create  a  test  data 
base.  The  first  consists  of  observations  of  range,  azimuth,  altitude,  and 
aircraft  identification  code  and  is  called  the  System  Demonstratior  !■  -'gram 
(SOP)  tape  (ref.  12).  Tne  other  source  is  derived  from  the  SDP  tipe  -nd  is 
valuable  for  editing  and  for  the  time-base  correction  which  is  described 
below.  This  latter  tape  is  referred  to  as  a  ZFLAG  tape  because  of  the  presence 
of  flag  wc.  .  which  are  used  to  qualify  each  target  report  (ref.  13). 

In  the  remaining  sections  of  this  Appendix,  a  brief  description  of  each 
of  the  main  data  preparation  steps  performed  on  the  SDP  data  will  be  given. 

These  steps  are  summarized  graphically  in  Figure  C-l. 

B .  Determination  o f  Accurate  Target  Report  Times 

The  time  associated  with  a  target  report  is  the  time  at  which  the  block 
itself  was  recorded  by  the  ATCRBS  target-to-tiack  correlation  task  and, 
therefore,  not  necessarily  the  desired  observation  time. 

In  order  to  rectify  this  situation,  and  generate  a  more  accurate  target 
tag  time  jequired  by  this  study,  a  procedure  was  developed  based  on  matching 
the  ’'average"  azimuth  vaiue  associated  with  a  target  report  on  the  ZFLAG  tape 
with  azimuths  of  individual  replies  on  the  SDP  tape.  Once  this  match  is 
achieved,  the  new,  more  accurate  time  is  obtained  by  linear  interpolation. 


--  gj 
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.  C-l.  Main  steps  in  data  preparation. 


Since  the  time  between  reply  blocks  is  10  ms,  the  calculated  report  time  is 
accurate  to  within  a  few  milliseconds  and  is  entirely  adequate  for  the  estimation 
of  a  time  bias  between  the  two  sensors. 

Shown  in  Figure  C-2  is  a  time  interval  during  which  the  aircraft  having 
a  code  of  46SS  was  executing  a  steep  ascent  maneuver.  After  the  time  has 
been  accurately  determined  for  both,  sensors,  a  time  offset  of  approximately 
1S00  ms  can  be  readily  discerned.  This  relative  time  offset  was  found  in  the 
data  from  other  aircraft  and  was  used  as  the  initial  best  estimate  of  the 
time  bias  TBI2. 

C .  Altitude  Correction 

The  altitude  recorded  on  both  the  Lexington,  Mass.,  and  Providence, 

R.I.,  data  tapes  is  not  corrected  for  barometric  pressure.  Fortunately, 
simultaneous  SDP  data  tapes  were  also  available  from  the  AKTS-IIJ  sensor  at 
Logan  International  Airport  in  Boston,  Mass. 

In  the  ARTS  system  the  altitude  recorded  in  the  target  report  block 
reflects  an  altitude  correction  not  made  to  the  altitude  appearing  in  the 
reply  block.  Thus,  the  difference  between  the  report  block  altitude  and  the 
reply  block  altitude  for  an  ARTS  target  is  the  adjustment  which  is  to  be 
applied  to  the  target  report  at  Lexington  and  Providence. 

It  was  found  that  an  altitude  correction  of  300  ft  had  to  be  added  to 
both  the  Lexington  and  Providence  data. 

D.  Conversion  to  Metric  Units 

In  this  final  phase  of  data  preparation,  the  target  report  records  are 
converted  from  units  used  on  the  SDP  tape  (ref.  12)  to  range  and  altitude  in 
meters  and  azimuth  in  radians.  Due  consideration  is  given  to  the  range  clock 
settings  at  each  sensor. 

Some  data  editing  is  also  performed.  For  example,  target  reports  with 
questionable  altitude  values  are  not  copied  into  the  final  test  data  base. 
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Fig.  C-2.  Example  of  time  base  offset  between  sensors 


APPENDIX  t) 

MATHEMATICAL  OVERVIEW 


A.  Int  roduct ion 

This  Appendix  includes  background  informat  ion  relevant  to  topics  included 
in  the  main  body  of  the  report.  The  level  of  coverage  is  not  intended  to  be 
rigorous  and  the  reader  interested  in  additional  detail  should  consult  references 
listed  in  the  bibliography. 

B.  Maximum  Likelihood  Bias  Estimator 

The  discrepancy  or  residual  in  apparent  aircraft  position  as  reported  by 

two  sensors  at  some  time  is  modelled  as  a  two-element  Gaussian  random 

vector  A . , 

1 


where  X  is  the  geographic  north-south  component  of  the  residual  and  Y  is  the 
east-west  component. 


The  probability  densitv  function  fui  A.  can  be  written  as: 
‘  1 


P(A.) 


1  /.T  _-l  .  , 
'2  (Ai  Ci  V 


|C.  |1/<:  2tt 


(DZ) 


T  - 1 

where  A.  is  the  transpose  of  A.  and  C.  is  the  inverse  of  the  covariance  matrix 

l  l  i 

of  the  residuals.  This  matrix  will  be  developed  below  on  the  assumption  that 
the  bias  values  are  known  and  that  the  residuals  arise  exclusively  from  measure¬ 
ment  error  in  both  sensors. 


For  N  independent  observations  of  the  aircraft,  the  joint  probability  density 
function  P(Aj,  A^,  .  .  •  A^)  is  given  by: 

N 

1  v  ,.T  _~1  .  , 

-T  )  (A.  c.  A. ) 

2  i  i  i 


P(A1(  A2,  .  .  .  y  = 


1 _ 

(271  )N  ^'=1|C.|1/2 


(W) 
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The  maximum  likelihood  choice  of  bias  parameters  is  based  on  the  idea 
that  the  underlying  system  should  be  the  one  which  makes  the  observed  set  of 
A^  much  more  likely  than  any  others.  Thus,  one  would  expect  the  best  estimate 
of  the  bias  parameters  to  be  that  which  maximizes  the  joint  probability  density 
function.  Because  of  the  negative  exponent,  this  amounts  to  minimizing  the 
function : 

N  . 

S  =  l  (A  j  C?1  A.)  (D4) 

i=l 

The  covariance  matrix  of  the  residuals  will  be  developed  by  first  con¬ 
sidering  the  measurement  error  covariance  matrix  from  one  of  the  two  radars.  Ke 
would  like  this  covariance  to  be  expressed  in  the  common  Cartesian  system  into 
which  all  local  measurements  are  transformed. 

From  the  geometry  of  Figure  B-l  it  follows  that  the  measurement  can  be 
written  in  Cartesian  coordinates  as: 


x  =  R  sin  9 
y  =  R  cos  0 

It  will  be  assumed  that  measurement  errors  in  R  and  6  are  each  zero  mean 
and  independent  and  have  a  covaraince  matrix. 


Now,  by  the  principle  of  the  propagation  of  covariances,  we  can  write. 


where : 
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Upon  performing  the  indicated  differentiations  and  matrix  multiplications,  the 
desired  measurement  covariance  matrix  for  one  sensor  can  be  written  as: 


C 

xy 


2  ■>  >  t  ? 

s :  n  0.  a~  *  R"  cos'"  6.  o„ 

l  R  l  10 


sin  6.  cos  Gi  (a2  -  R2  a2) 


sin  6.  cos  0.  (o“  -  R."  o7)  cos"  6.  of  +  P."  sin"  0.  of 
l  i  R  10  lR  l  10 


Residuals  in  the  x  and  y  directions  can  be  simply  expressed  in  terms  of 
the  latitude  and  longitude  (<£,  X)  of  the  aircraft's  position.  Specifically, 


x  =  E  cos  <f>2  (X  -  Xj ) 

y  =  E  (<t>2  -  Oj) 

where  L  is  the  mean  radius  of  the  referenced  earth  ellipsoid  and  the  subscripts 
on  0  and  X  identify  the  reporting  sensor.  The  factor  cos  in  the  expression  for 
x  was  arbitrarily  selected  as  a  scaling  factor;  cos  4,  could  have  been  used 
instead . 


Recall  that  A.  is  calculated  from  the  difference  between  two  observations, 
i 

By  taking  these  to  be  independent,  the  covariance  of  the  residuals  is  then  just  the 
sum  of  the  covariances  of  the  measurement  errors  from  the  two  sensors.  If  this 
sum  matrix  is  written  as: 


..sensor  1  ..sensor  2 
C.  =  C  +  C 

i  xy  xy 


or 


r  = 
"i 


bll  b12 


b12  B22 


then  the  function  S  to  be  minimized  can  be  expressed  as: 


S  = 


1 


2  1  22 
(bH  b22  -  bj2) 


(b->,  x2  -  2  b12  xy  +  bn  y2) 
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It  can  be  shown  by  substitution  that  S  is  the  sum  of  the  squares  of  the 
scaled  projections  of  the  vector  A  on  the  major  and  minor  axes  of  the  error 
ellipse.  The  scale  factors  are  the  error  ellipse  axes  themselves.  Thus,  the 
length  of  a  difference  vector  alone  is  not  as  important  in  determining  the  biases 
as  is  the  length  and  orientation  with  respect  to  the  combined  error  ellipst  at 
the  aircraft  position. 

Presented  in  Figure  D-l  is  a  geometrical  interpretation  of  the  term 
C. 1  at  the  i*"*1  "superpoint"  sample  time.  As  the  aircraft  changes  position 
with  respect  to  the  two  sensors,  the  effective  weighting  of  a  residual  vector 
of  a  given  lengch  varies  accordingly. 

C.  Variance  of  the  Bias  F-stimate 


Analysis  of  the  error  associated  with  the  bias  estimates  considered  in 
this  report  would  be  quite  complex  owing  to  the  non- linear  character  of  the 
problem.  It  is,  however,  instructive  to  consider  instead  a  simpler  problem: 

Given  a  set  of  observations  x  ,  y.  i  =  1,N,  and  an  assumed  function 

l  i 


y(x)  =  l  a.  Q  (x) 

j  =  0  J  J 

Find  the  parameters  a^  such  that  the  quantity 


{>i  -  y  (xL)}2 

c. 

1 


is  minimized. 


(x)  is  typically  a  polynomial  in  x  while  is  the  sample  variance  of 

the  y .  '  s . 

'  i 

Required  in  the  development  of  the  variance  of  the  parameter  estimates  is 
the  K  by  K  symmetric  matrix  H  whose  elements  re ^  are: 

i  a2  x2 

jll  2  3a.  3a. 

1  X/ 
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.  D-l.  Geometrical  interpretation  of  A  C  A. 
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This  matrix  is  referred  to  as  the  Hessian  matrix  and,  more  descriptively, 

2 

as  the  curvature  matrix  because  of  its  relationship  to  the  curvature  of  X  in 
parameter  space. 

It  can  be  shown  that 


*» '  Jo { j,  l7*  yi  ^  <x*: 


where  is  an  element  of  the  matrix  E  where 


E  =  H 


The  uncertainty  in  the  parameters  can  be  estimated  from: 


-  N  -  3a. 

°a.  =  J  ai 

J  i=l  }x 


The  derivative  in  the  summation  can  be  expressed  as: 

3a.  K  e.j 

*  Jo 

1 

Substituting  and  rearranging,  the  variance  of  the  parameter  estimate 


becomes : 


a  =  e . . 
a.  j  j 

Because  of  this  relationship,  the  matrix  E  =  (e^)  is  called  the  error  matrix. 

D.  Combining  Biases  from  Two  or  More  Aircraft 

Having  obtained  an  estimate  of  a  bias  a  from  the  m^  aircraft  and  its 
2  ^ 

associated  variance  o  ,  it  will  now  be  shown  how  estimates  from  M  different 
aircraft  should  be  combined  to  obtain  the  best  possible  composite  estimate  a. 
Assuming  Gaussian  probability  again,  we  maximize  the  probability  of  getting  the 
observed  M  estimates  by  minimizing  the  associated  probability  density  function. 
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The  exponent  of  this  function  can  be  shown  to  be: 

i  »  <*b-»a 

' 7  Ji  “v7- 

am 

and  thus  the  task  at  hand  is  to  minimize  by  taking  3/3am  and  setting  equal  to 
zero. 


1  3 


M  (a  -  S)‘ 


2  JT 

m  m=l 


simplifying 


M  (a  -  S) 
v  m 

2.  9 


solving  for  3: 

M  a  M  . 

a  *  I  -f-/  I  -4- 

m=l  o  m=l  a 

a  a 

m  m 

Thus,  the  composite  bias  estimate  is  a  weighted  mean  and  its  variance  is 
given  by 


n  i 

l  (—Vo 


m=l  a 


E.  Mean  and  Variance  for  S 


Given  the  statistic  S 


,  N  x  . 

s  s  I  y  __L_ 
v  ■  S  „  2 
1=1  a 

x 


where  x.  is  assumed  to  be  a  zero  mean  Gaussian  variable  and  x.,  x,  for  j  f  k, 
i  1  "  2 

are  independent.  Also  assumed  is  that  the  variances  are  equal,  i.e.,  c r  = 

d2  for  all  K.  It  is  well  known  that  S  has  a  chi-square  distribution  wii:h  v 


degrees  of  freedom. 
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Expressions  for  the  mean  and  variance  of  S  will  now  be  developed. 


Mean 


N  x? 


E  [S]  ■  E  [±  l  — y- ] 
v  i=l  cf 


i  N  ? 

“4"  E  I  l  xi  1 

v  ox  i-1 


v  a2 


N 

\) 


Variance 


Var(S)  «  E  [(S  -  E[S]H 

N  x2 


*  E  [(i  l  -i-  -  E  [S])2] 


i=l 


Substituting  for  the  mean  and  expanding: 

N  x?  ,  N 


N  x 2. 


„  _  r,l  r  *1  w1  V  %  2  r  V  L  A21 

*ar(S)  E  [(  I  ?  V  2^  v  ^  ^  2^  v  V  ^ 

j"1  <5  K,=1  °x  1=1  % 

By  carrying  through  the  expectancy  operation  for  each  term, 


N  x2 


N  x 


n  rr1  V  r1  V  %,  P  r2N  V  Kl  *  PrN  i 

s  E  It:  1  -rJ  (u  i  -t)]  *E  hi  -j]  +  E[— jj 


'V  .*•,  2 

j=l  a 
J  x 


Vv  **.  2J 
K-1  ox 


The  middle  term  becomes: 


2  tNc?x  ) 

■  v  c~r~  J 


N 

v 


2N 


and,  thus,  the  second  and  last  terms  combine  to: 


75 


-N2 

What  remains  then  is  to  reduce  the  first  torm;  the  expectancy  of  the  product  of 
two  summations. 

Factoring  out  the  from  the  denominator  in  the  summation,  the  first  term 
can  be  written  as: 

N  „  N 


-r~r  E  (<  l  1 1  x2)) 

v  or  j=i  3  k=i 


The  product  of  the  summations  includes  w  terms  of  which  N  are  of  the 


form: 


x4 

3 


3  -  l.N 


2  2 
Xi  Xj 


and  it  can  be  shown  by  Gaussian  moment  factoring  that 

E[x4]  .  3  a4 

2 

The  remaining  N  -  N  terms  of  the  product  are  of  the  form 

i  f  3 

and  since  the  x^'s  are  independent  and  have  equal  variances  and  zero  mean 

E  [x2  x2)  =  E  [x2]  E  [x2]  '  oj;  =  04 

Thus,  the  term  above  involving  the  product  of  the  summations  becomes: 

1  „  x,  J  .  ,.,2  _  N2  ♦  2  N 

— 5 — 2 —  (3  N  o  +  (N  -  N)  o  )  =  - ~ - 

v  a  x  x 

x 

and,  finally,  the  variance  of  S  can  be  obtained  by  combining  results  and 


Var(S)  = 


N  +  2  N 


N 


2  N 
.  2 
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For  the  problem  considered  in  this  report,  v  is  the  number  of  degrees  of 
freedom 

V  *  N  -  L  -  1 

where  L  is  the  number  of  bias  components  present  and  N  is  twice  the  number  of 
"superpoints"  employed  in  the  estimation  (one  for  each  of  the  x  and  y  residuals) . 

For  N  large,  the  mean  of  S  approaches  1  while  the  variance  approaches 

2/N. 
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Today’*  tactical  military  air  surveillance  radars  generally  operate  in  a  stand-alone  configuration.  The  many  performance 
improvements  that  result  when  data  from  multiple  radars  of  this  type  are  merged  have  made  such  netted  operations  an  “tractive 
goal  for  many  years.  A  major  obstacle  to  achieving  this  goal  has  traditionally  been  the  difficulty  associated  with  the  registration 
of  multisensor  data,  the  expression  of  the  data  in  a  common  coordinate  system  free  from  errors  due  to  site  uncertainty,  antenna 
orientation,  and  Improper  alignment. 

This  report  presents  the  resolts  of  a  modest  effort  to  develop  a  self-registration  procedure  by  which  multiple  radar  sensors 
operating  in  consort  each  calculate  the  errors  in  tlielr  data  by  comparing  it  with  data  from  the  remainder  of  the  system  and  then 
uses  the  Information  to  upgrade  performance.  The  technique  has  been  tested  with  experimental  data  and  appears  quite  capable  of 
improving  system  performance,  measured  in  terms  of  residual  inter-site  bias  errors,  by  almost  a  factor  of  one  hundred. 
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